Decoupled and unidirectional asymptotic models for the 
propagation of internal waves 

Vincent Duchene* 
September 3, 2012 

Abstract 

We study the relevance of various scalar equations, such as inviscid Burgers', Korteweg-de 
Vries (KdV), extended KdV, and higher order equations, as asymptotic models for the propa- 
gation of internal waves in a two-fluid system. These scalar evolution equations may be justified 
in two ways. The first method consists in approximating the flow by two uncoupled, counter- 
propagating waves, each one satisfying such an equation. One also recovers these equations 
when focusing on a given direction of propagation, and seeking unidirectional approximate 
solutions. This second justification is more restrictive as for the admissible initial data, but 
yields greater accuracy. Additionally, we present several new coupled asymptotic models: a 
Green-Naghdi type model, its simplified version in the so-called Camassa-Holm regime, and a 
weakly decoupled model. All of the models are rigorously justified in the sense of consistency. 

1 Introduction 

In this paper, we study asymptotic models for the propagation of internal waves in a two-fluid 
system, which consists in two layers of immiscible, homogeneous, ideal, incompressible and irro- 
tational fluids under the only influence of gravity. We assume that there is no topography (the 
bottom is flat) and that the surface is confined by a flat rigid lid, although the case of a free surface 
could be handled with our method. Finally, as we will focus on scalar models, we restrict ourselves 
to the case of horizontal dimension one. The precise equations of such a system are presented in 
Appendix |Aj and we refer to the survey article [31] and references therein for a good overview of 
the relevance and significance of such system in oceanography. 

The mathematical theory for the governing equations of our system is extremely challenging 
(see [43] in particular), which has led to extensive works concerning simplified, asymptotic models, 
aiming at capturing the main characteristics of the flow with much simpler equations, provided the 
size of given parameters (typically measuring the amplitude of the deformation at the interface, 
and/or shallowness of the two layers of fluid) is small. Among other works, we would like to 
highlight 8J, where many asymptotic models are presented and justified, in a wide range of such 
regimes. In each case, the resulting model consists in two relatively simple evolution equations 
coupling the shear velocity and the deformation at the interface. 

However, in this work, we are mainly interested in decoupled scalar models, which can be used to 
describe the unidirectional propagation of gravity waves. The derivation and study of such models 
have a very rich and ancient history, starting with the work of Boussinesq |llj and Korteweg-de 
Vries [41] which introduced the famous Korteweg-de Vries equation 

(KdV) d t u + cd x u + aud x u + vd x u = 0, 

in the context of propagation of surface gravity waves, above one layer of homogeneous fluid (which 
we refer as "water-wave problem"). However, let us note that the complete rigorous justification of 
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such model is recent: [571 [S5J [7J. Such a justification is to be understood in the following sense: in 
the long wave regim^ 

e ~ /z < 1, 

the flow can be accurately approximated as a decomposition into two counter-propagating waves, 
each component satisfying a KdV equation (or more generally, with the same order of accuracy, a 
Benjamin- Bona-Mahony equation [51 154]). 

(BBM) d t u + cd x u + aud x u + v x d x u — u t d 2 d t u = 0, 

The situation is similar in the case of two layers of immiscible fluids, with a rigid lid (as well as in 
the case of two layers with a free surface, except the flow is then decomposed into four propagating 
waves), and correspondent results are presented in [23| . 

Of course the coefficients (c, a, v, etc.) depend on the situation, and a striking difference be- 
tween the water-wave case and the case of internal waves is that in the latter case, there exists a 
critical ratio for the two layers of fluid (depending on the ratio of the mass densities) for which 
the nonlinearity coefficient a vanishes. In that case, heuristic arguments support the inclusion of 
the next order (cubic) nonlinearity, which yields the modified KdV equation (see [52| and references 
therein): 

(mKdV) d t u + cd x u + c<2U 2 d x u + v x d x u — v t d 2 d t u = 0. 

For the cubic nonlinearity to be of the same order of magnitude as the dispersive terms, this urges to 
consider the so-called (refering to |19) for this non-standard denomination) Camassa-Holm regime: 

e 2 ~ fi < 1. 

Note that this regime is interesting by itself as it allows larger amplitude waves than the long wave 
regime, and in particular yields models developing finite time breaking wave singularities [181 119) . 
Thus we ask: Can we extend the results of the long wave regime to the Camassa-Holm regime ? 

As we shall see, the answer is not straightforward, as several evolution equations are in compe- 
tition, with an accuracy depending on the situation, and in particular the criticality of the depth 
ratio, and localization in space of the initial data. In addition to the equations already presented 
above, we consider the extended KdV (or Gardner) equation 

(eKdV) dfU + cd x u + aiud x u + ct2U 2 d x ii + v x d x u — i> t d 2 dtu = 0, 

as presented in [36l EH 06] , studied in [29l [5j2 [55] (among other works) , and tested against ex- 
periments in [40]. l32[ [45] . More generally, we work with the higher order model presented in the 
water-wave setting in |19j (which we therefore call Constantin-Lannes) 

dtu + cd x u + a\ud x u + a^u d x u + a^u d x u + v x d 3 x u — u t d 2 dtu 
(CL) + d x (Kiud x u + K2(d x u) 2 ) = 0, 

which is related, though slightly different from the Camassa-Holm equation obtained for example 
in [T31 [35j or higher order models (see various such models in [SU [551 123 E01 \EH I2H HO] , concerning 
the water wave or internal wave problem). However, each of these works is limited to the restrictive 
assumption that only one direction of propagation is non-zero. 

We propose to study in detail the justification of the decomposition of the flow (although the 
unidirectional case is also treated) as presented above, in the case of internal waves. Our results 
are expressed in a very general setting: assuming only 

fi <C 1 and e <C 1, 

and expressing the accuracy of the different models as a function of these parameters. Although this 
complicates the expression of our results, and renders the proofs fairly technical, such choice allows 
to cover general regimes, including the long-wave as well as the Camassa-Holm regimes, described 
above. We recover in particular the relevance of the KdV approximation in the long wave regime. 
In the Camassa-Holm regime, the situation will depend on the criticality of the depth ratio, the 
localization in space of the initial data, as well as the time-scale which is considered. 



1 in the following, fi is the shallowness parameter, and e is the nonlinearity parameter, as defined precisely (in the 
bi-fluidic case, see [7] for the water-wave case) in Appendix [A| 
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Outline of the paper, and overview of the results. For the sake of readability, the precise 
presentation and justification of the governing (full Euler) equations of our system is postponed to 
Appendix |A| Using the shallowness assumption (/i <C 1), we then introduce the so-called Green- 



Naghdi model (1.1), written out below. This system is related to the one introduced by Choi and 
Camassa in |16] ~ We present its first rigorous justification, in the sense of consistency: roughly 
speaking, any sufficiently smooth and uniformly bounded solution of the full Euler system satisfies 



the (A. 19) model, up to a residual of size 0{ji ). Our decoupled models will be constructed from, 
and compared with the Green-Naghdi model. 

In section |2.1[ we present a formal argument which allows to derive the decoupled models at 
stake, that we define precisely in Definition |2.3| and Definition 1 2. 4 1 

The main results of the paper are presented in Section [272] We first show in Proposition [275] that 



the equation (CL), and therefore any simplified equation, is well-posed and uniformly controlled 
in Sobolev norms over long times, for sufficiently smooth initial data. We also prove that the 
localization in space of the initial data (that we express by the fact that the function lives in 
weighted Sobolev norms) propagates over long times. 

We express our main result in Proposition |2.7| and Corollary |2.10[ stating that for any suffi- 
ciently smooth initial data, an approximate solution of the Green-Naghdi equation (in the sense of 



consistency) consists into two decoupled waves satisfying ( CL ) (or a lower order equation) , plus an 
explicit, small coupling term. The precision in the sense of consistency and estimates on the size of 
the coupling term depend on 

i. the evolution equation considered (CL,eKdV,KdV, etc.); 

ii. the size of the parameters e, /i as well as 5 2 — 7 (critical ratio); 

iii. the localization in space of the initial data; 

iv. the time-scale considered. 

As announced above, our result is as extensive as possible, and a discussion on several impor- 
tant cases (long-wave regime, Camassa-Holm regime with critical or non-critical ratio) follows in 
Section [2.3[ with several numerical simulations to support our conclusions. 

As intermediary steps for the justification of decoupled asymptotic models, we introduce two 



coupled asymptotic models, which are interesting by themselves. The first one, ( 2.3a )-( 2.3b) is a 



simplified version of the Green-Naghdi equation, with same precision provided that we are in the 



Camassa-Holm regime (precisely defined in Remark 2.1), as stated in Proposition 2.2 We also 
introduce weakly coupled models, which add first order corrections of the fully decoupled models, 
in the spirit of |62) . These corrections are explicitly obtained using the solutions of the decoupled 
models; see Remark 12.81 



In Section [3] we turn to the case of a unidirectional flow. More precisely, and in the spirit 



of |35 | ll9 j . we show in Proposition 3.1 that one can construct very precise unidirectional approximate 



solutions: the interface deformation satisfies a Constantin-Lannes equation (CL), and the shear 
velocity is given as a function of the interface deformation. As we follow a specific direction of 
propagation, the coupling terms which were the main contributors for the error in the previous 
section vanish, and the accuracy of the approximation is therefore considerably improved. We then 
numerically investigate if the quite restricting condition on the initial data arises naturally, that is 
if after some time, the flow generated by any initial perturbation will eventually decompose into 
two almost purely unidirectional waves. 



Finally, Appendix [B] contains the proof of Proposition 2.5 and Appendix [C] the proof of Propo- 
sition 12.71 



The Green-Naghdi coupled model. Let us display here the Green-Naghdi model (A.19), as 
derived in Appendix [X] The system has two unknowns: C, representing the deformation at the 
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interface, and v the shear layer-mean velocity, as defined precisely in (A. 18) 
h±h 2 



dtC + dj^ 



(1-1) 



jh 2 



v = 0, 



d t (v + vQ[hiMv) + (l + t)d x ( + l d *( { h 1+ ^h 2 )^ 2 ) = M^x(^i,/i 2 ]w), 

where 7 and S are, respectively, the mass density and depth ratio (see Appendix [A} ; hi = 1 — eC 
and h 2 = k + eC are the depth of, respectively, the upper and the lower layer, and where we define 
the following operators: 



Q[hi,h 2 ]V 
n[h u h 2 ]V 



1 



V 



3 hi + 7/12 \ h 2 



-djh^d^- 
2 V h x 



h 2 V_ 
jh 2 

hiV s\ h 



lh 2 



h 2 V 
lh 2 



This model is justified the consistency of the solutions of the full Euler system towards the Green- 
Naghdi models, as expressed precisely in Proposition |A.5| 

As it has been said above, our models will be justified with consistency results, with respect to 



the Green-Naghdi model (1.1 1. Let us now define precisely what we denote by consistency in the 
core of this paper. 

Definition 1.1 (Consistency). Let (£ p ,v p )p^-p be a family of pair of functions, uniformly bounded 
in I/°°([0, T); H s+S ) (s > to be determined), depending on parameters 

(1.2) peV = {(e,M,<5,7), 0< fj, < /w, < ), < 7 < 1 }. 



We say that (C p iV v ) is consistent with Green-Naghdi system (1.1 1 at precision 0(e p ), of order s 
and on [0,T), if(C v ,v' 1 ') satisfies, fore' 1 ' sufficiently small, 



hi + jh 2 



d t (vP + vQ{hi,h 2 y) + h + s)d x c p + ^( ^^J ^l 2 ) ~ ^d x (n[hi,h 2 y) = r 2 , 

where hi = 1 — e£ p and h 2 — I + e£ p , and with 



I (j*l 1 ^2 ) J I ^co qq j*y jjs )2 — C ^A^max? ^maX7 f /^min ; ^maxi \\C^ 



l£°°([0,T);ff 8 + 3 )' 



li~([0,T);ff»+ 5 ) 



Here, and in the following, we denote by C(Ai, A2, . . .) any positive constant, depending on the 
parameters Ai,A2,..., and whose dependence on Xj is assumed to be nondecreasing. Moreover, 
for < T < 00 and f(t,x) a function defined on [0,T] x E, we write / € L°°([0,T];H s ) if / is 
uniformly (with respect to t G [0, T]) bounded in H s — H S (M.) the L 2 -based Sobolev space. Its 



norm is denoted 



Il°°([o t)-H"V as ^ e Sobolev norms are denoted with simple bars: | ■ \ R3 



2 Decomposition of the flow 



In this section, our aim is to obtain approximate solutions of ( 1.1 1, through a decomposition of the 



flow into two independent waves, each one satisfying a scalar evolution equation. Our aim is dual. 
First, we want to investigate which scalar equation each of these waves should to satisfy, in order to 
be as accurate as possible. Then, we want to estimate the size of the error we commit by neglecting 
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the coupling between the two components. Section [3] is focused on the case where coupling terms 
vanish, as only one of the components is non-zero. 

The strategy presented here has been used by [7j, where the authors present a rigorous justifi- 
cation of the KdV equation for the (one fluid) water wave problem in the long wave regime, and 
in [23] in the bi-fluidic case. Our work extends the latter results to more general regimes and scalar 
equations. 

We first give a fairly simple formal approach in section 2.1 which allows to heuristically construct 
the decoupled equations at stake, and consequently the approximate solutions. These approximate 
solutions are precisely defined in Definitions |2.3| and |2.4| In section |2.2( we give the rigorous 
justification of such approximations; see Proposition |2.7| and Corollary |2. 10 Section [273] contains a 
discussion on our result, considering various different regimes and decoupled models, and supported 
with numerical simulations. 



2.1 Formal approach 

The main idea of the decomposition is that at first order (that is, setting e = fj, = 0), our system of 
equations (1.1) is simply a linear wave equation 

(2.1) d t U + A Q d x U = 0(e,/i), with A =(° +s ^ 

and U = (C, v) T ■ It is straightforward to check that Aq has two distinct eigenvalues, therefore, we 
can find a basis of K 2 such that (2.1) reduces to two decoupled equations. More precisely, there 
exists 



P = 



1 



1 



j + S -(7 + 5) 



• P- 1 = 



f A 



tff such that P^A Q P = 
1 i+sJ 



1 

-1 



Now, let us define (ui,u r ) = P^U = §(( + C - ^)- Multiplying (2TJ) by P- 1 (on the left), 
and keeping only first order terms, yields 



(2.2) 



di 



1 

-1 



0(e,n). 



As a conclusion, the solution U — (£, v) T may be decomposed in the following way: 
U?* (v+(x-t) + v-(x + t) , (i + 6)(v + (x-t) + v-(x + t))^ . 



Now, let us take into account the higher order terms in ( 1.1 ). One will make use of the following 
straightforward expansions: 



hih 2 



1 



5 2 



7 



hi + jh,2 7 + 5 
hi 2 - 7/12 2 I 



2 7*0? + I) 2 ,2 



{hi+jh 2 ) 2 
Q[hx,h 2 ]v = 

H[hi,h2;]v = 

1-7 



(7 + *) 



7 _ 2c 7^ + l) 2 c _ 3e2 7* 2 (* + l) 2 (l-7) 



1 + 7* 02- 



( 7 + *) 3 
7 + * 



7* 2 (* + l) 2 (l- 7 ) 
(7 + *) 4 



C 3 + 0(e 4 ), 



afyd x v? + \vdlv\ + 0(e) 



with 



and (3 = 



(l + 7*)(* 2 -7) 
*(7 + <5) 3 ' 



Using the decomposition as above: (ui,u r ) = P 1 (( } v) T = ^((+ , ^g)> and withdrawing 
every term of size 0(fie 2 , e 4 ), one can check that the Green-Naghdi system ( |1.1[ ) becomes the 
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following coupled system in terms of ui (left-going wave) and u r (right-going wave): 

(2.3a) d t ui + d x ui + fl^(u u u r ) = 0, 

(2.3b) d t u r - d x u r + f^(u h u r ) = 0, 



where / ; £,AI and / ; e,At are defined below: 



J i 



e ^-9 x ((ui + -u r )(ut - U r )) + e 2 ^-d x {{ui - u r )ui(ui + u r )) 



,3^3. 

4 



+ e°— d x ((ui - ~u r )(ui - u r )(ui + u r ) 2 ) 



[ivd 2 x d t {ui - u r ) + [ieK 3 d x (^(ui - u r )(d 2 ui - d 2 x u r ) + ^(d x ui - d x u r ) 2 ^j 



fiedt ni{uidlui - u r dlu r ) + K 2 {u r dlui - uid 2 u r ) + (ki + K ^-){{d x ui) 2 - (d x u r f 



fr 



-e^-d x ((^-ui + u r )(u r - ui)) - e 2 ^-d x ((u r - Ui)u r (m + u r )) 



2 vv 3 
e 3 ^-d x ({u r 



~ : ui)(u r - ui)(ui + u r y 



[ivd 2 x d t {u r - ui) - neK 3 d x (^(u r - ui){c^u r - d 2 u t ) + ^(d x u r - d x ui) 2 ^j 



fiedt ni(u r dlu r - uidlui) + n 2 (uidlu r - u r dlui) + (ki + ^-)({d x u r 



(d x ui) 2 



where we used the notation 

a , - 3^ 
(2.4) 1 2 i +s ' 



_ q 7 <5(,5+1) 2 
" 2 — (7+<5) 2 ' 
_ (l +7 g)(g 2 - 7 ) 
— 3<5( 7 +<5) 2 ' 



_ ,5 2 (,S+l) 2 7 (l- 7 ) 

a 3 - -5 ( 7 +yp , 

_ (1 — y) _ 7 -l 

K2 - 3( 7 +<5) ' K 3 - 2(7+5 ) ■ 



__ 1 i+ 7 a 

— 6 <5(<5+ 7 ) ■ 



Remark 2.1. One could use higher order expansions with respect to the parameter e, which would 
lead to decoupled models with formally higher accuracy. However, let us note that our results (see 
discussion in section 2.3) show that the main error of such decoupled models comes from neglecting 
coupling terms which arise at low order, so that including higher order terms in the evolution 
equation is unlikely to produce substantial improvement. 

Recall that the Green-Naghdi system is consistent with the full Euler system at precision 0(^l 2 ); 
Proposition A. 5 Therefore, choosing to keep only terms of order C(/xe 2 ,e ) in (2.3a)-(2.3b) is 
especially appropriate for parameters in the so-called Camassa-Holm regime: 

(2.5) V C h = {(e,n,5,i), 0< fi < /x max , < e < Myfjx, 6e (5 

which is natural to consider, as discussed in the introduction. 



mini '-'max 



), 0< 7 <!}, 



Proposition 2.2 (Consistency of (2.3a )-( 2.3b)). Let (itf,uf.) be strong solutions of (2.3a|-(2.3b|, 
depending on sets of parameters (e, fj,, 8, 7) — p G V, as defined in (1.2 1. We assume that for 
s > s > 1/2, u\,v? r £ W 1 {[0 1 T);H s+3 (R)), uniformly in p, where Wlfc, T); H S+1 (R)) denotes 
the space of the functions f(t,x) such that 



l/l 



/ 



Additionally, we assume that there exists h > such that 



I ft/I 



L°°([0,T);/f 5 



< OO. 



h{ = 1- e(u£ +uf.) >h>0, 



h p 2 = - + e{u v l + u?) >h>0. 



Then (C p ,w p ) = K +u», (7 + S)(u^ — uf.)) is consistent with Green-Naghdi equations ( |l.l| of 
order s on [0,T), at precision 0([j,e 2 + e 4 ) (in the sense of Definition 1 1 . ip . 

Proof. Proposition |2.2| is straightforward once one obtains a rigorous statement of the expansions 
°f hl+^h" ' ^(feP +-yhly ' Q-VAi^l]^ 1 T^-VAi hl)^ , as stated above, the residual being estimated 
in W 1 ([{),T);H S+1 ) norm. These expansions are easily checked, provided the assumptions of the 
proposition (u^uf. £ W Al ([0, T); £P +3 (R)), h\,h\ > h > 0) are satisfied, and using uniformly 
continuous (for s > s > 1/2) Sobolev embedding H s <-t L°° . □ 
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The uncoupled approximation consists simply in neglecting all coupling terms in (2.3a I— ( 2.3b I , 
that is replacing f^(ui,u r ) by /f' M (uj,0), and f^ l {u h u r ) by f^(0,u r ). This yields 



(2.6a) 



(2.6b) 



dm + d x ui + eaiuid x ui + e 2 a 2 ufd x ui + e 3 a 3 u 3 d x ui - fii/d x d t ui 
-/j,edt(niuid%.ui + («i + l/2n 2 ){d x ui) 2 ) + [ien 3 d x (^uidlui + ^(d x ui) 2 



dtu r — d x u r — ea\u r d x u r 



2 2a 

e a 2 u r o x u r 
\2\ 



e 3 a 3 u r 3 d x u r — iivd 2 dm 



-Lied t (Kiu r dlu r + (ki + l/2n 2 ){d x u r ) 2 ) - fieK 3 d x {^u r dlu r + ]^{d x u r ) 2 ) 



= 0, 



0. 



and ( 2.6a H 2.6b I are the decoupled equations we consider; see Definition 2.3 below. 



Let us now reckon that one can deduce from ( 2.6a H 2.6b I a large family of formally equivalent 
models, with different parameters, following the techniques used for example in [6j El HH [23] , and 
that we discuss below. 

• The BBM trick (from Benjamin-Bona-Mahony [5 ). Keeping only the first order terms in 
equations ( 2.6a[ )-(2.6b|, one has the simple transport equations 



dm + d x ui = 0(n, e), and d t u r 



d~ 



It follows that one can replace time derivatives in higher order terms (of order O(pe)) by 
spatial derivatives (up to a sign), and both equations have formally the same order of accuracy. 
Throughout this paper, we consider only equations with spatial derivatives in O(fie) terms 
(if they exist), to simplify the rigorous approach of Section 2.2 In particular, (2.6a)-(2.6b) 
becomes 



= 0, 



dm + d x ui + eamdxUi + e 2 a 2 ufd x ui + e 3 a 3 u 3 d x ui - [ivd 2 x dm 
(2.7a) +Hed x ((K! + 1/3^)^6^ + (« a + 1/2k 2 + l/2 K3 )(d x inf) 

dm — d x u r — ea\u r d x u r — e 2 a 2 u 2 d x u r — e 3 a 3 u r 3 d x u r — [ii>d 2 dtu r 
(2.7b) -fied x ((Ki + l/3n 3 )u r d*u r + (k x + 1/2k 2 + l/2n 3 )(d x u r ) 2 ) = 0. 



Following a similar idea, we make use of the low order identity obtained from ( 2.7a |— ( 2.7b I : 

dm + d x u t + eamid x ui = 0((i,e 2 ), 
so that one has, for any 9 £ K, 

dm = Qdm + {9 - l)(d x ui + ectiuid x ui) + 0{p,,e 2 ). 



Plugging back into (2.7a) and withdrawing C(/i 2 ,/xe 2 ) terms yields 

dm + d x ui + ea x uid x ui + e 2 a 2 ufd x ui + e 3 a 3 u 3 d x ui - iiv e t d 2 x dm 
(2.8a) + f^ x d 3 x ui + ^d x (4uid 2 x ui + 4{d x ui) 2 ) = 0, 

dm — d x u r — ta\U r d x u r — e 2 a 2 u 2 d x u r — e 3 a 3 u r 3 d x u r — ^,v^d 2 dm 
(2.8b) - \i x vd 3 x u r - ^d x {K x u r dlu r + K%(d x u r ) 2 ) = 0, 



where we have defined, after parameters (2.4 1, 

(2.9) vl = Ou, vi={l-6)v, k\ = Ki + y + (l-0)arii/, 4 = ^i + y + y + {l-9) ai v. 

• Near identity changes of variables. We used system ( |1.1[ ) as our reference system, and therefore 
the unknowns we consider are zl) , where v — u 2 — jui is the shear layer- mean velocity. 
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However, other natural variables may also be used, such as the the shear velocity at the 
interface (leading to system (A. 17)) 

va = d x ((^2 - 701 ) |*=ec) = d x (j) 2 (x,r 2 (x, 0)) - j(f>i(x, n(x, 0))), 

(where we use the change of coordinate flattening the fluid domains: ri(x, z) = zhi(x)+e£(x); 
see [12] ), or more generally, using the horizontal derivative of the potential at specific heights 
( Zl ,z 2 ) E [0,1) X (-1,0]: 



d x (4> 2 {x,r 2 (x,z 2 )) - 7^1 (x, r 2 (x, z\)) 



Using the expansion of the velocity potentials and layer-mean velocities, as obtained in |22j . 
as well as the identity h\Ui + h 2 u 2 — (obtained through the rigid lid assumption), yields 
the following approximation: 



= v + fj,X zl ' Z2 d 2 V + fieT Zl ' Z2 (C,v) + 0(n 2 lt j,e 2 ), 

with \ Zl ' Z2 = — — — ® Z2 — "^ff 1 1 — — -, and f Zl ' Z2 a bi-linear differential oper- 

65(7 + S) 

ator whose precise formula do not play a significant role in our work as we shall discuss below. 



Following this idea, we consider 



ui + /j,\d 2 ui 



.A _ 



= u r — /iXd x u r 



If (ui,u r ) satisfies (2.8a)-(2.8b), then (u x ,u x ) satisfies the following equations, up to terms 
of order 0(fi 2 , ^e 2 ): 

d t u x + d x u x + eaiu x d x u x + e 2 a 2 u x2 ' d x u x + e 3 'a^uf 3 ' d x u x — fi^' X d x d t u x 



(2.10a) 
d t u x 
(2.10b) 



+ ^ e x x d 3 x u x + (xed x (4MdW + K e 2 (d x ut) 2 ) =0, 

n A A n A 2 A 2 n A 3 A 3 o A 6, A a2 o A 

o x u r — tct\u r o x u r — e a 2 u r o x u r — e a^u r o x u r — fxv t ' a x o t u r 

- lu4- x B%v£ - ^ed x {^ X u x d 2 x u x + K 6 2 {d x u x ) 2 ) =0, 



where we have defined, after parameters (2.9|, 

0,A _ e , \ e, a _ 8 \ 0.A _ e.A , \ 



(2.11) 



Note that in order to fit as much as possible with variables ((,v zl ' Z2 ), one could have used 
more complex change of variables, such as 

uf = ui + ^i\d' 2 ui + ^ie(X 2 uidlui + \ 3 d x (u 2 )), 
u x . = u r — /iXd 2 u r — /ie(X 2 u r d 2 u r + X^d x {u r 2 )) , 

with A = 2(j+s) ' anc ^ ^2,^3 obtained through 7~ z i> z 2. It turns out u x and u x would then 
satisfy the same equation as above: the new parameters do not depend on A2 and A3, as their 
contribution is of order C(/xe 2 , /i 2 e), after using BBM trick to suppress higher order derivatives 
with respect to time. We thus do not consider such changes of variable. 



Ultimately, one obtains the following family of approximations: 
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Definition 2.3 (Constantin-Lannes approximation). Let be given scalar functions, and set 

parameters (e,/i,7,5) G V (see (1.2)), (A, (9) € R 2 . The Constantin-Lannes approximation (CL) is 

U C L = (v + (t,x-t)+v-(t,x + t),('j + 5)(v + (t,x-t)-v-(t,x + t)) > ), 

where v± \t=a = ^(C° ± ^+s) \t=o an d v± = (1 ± ii\d 2 )~ 1 v± with v± satisfying 
(2.12) d t vi ± e ai v±d x vl ± e 2 ^^) 2 ^ ± eS^'^) 3 ^ 

± ^ A dM - iiv!' x %dtvi ± ^(^' A ^a 2 ^ + ^(a^i) 2 ) = o, 

wi£/i parameters defined as 

3 5 2 - 7 7<5(J + 1) 2 5 2 (5 + 1) 2 7(1 _ 7) 

= „ —r, ct2 = -i>—, — , rN0 ) "3 = -5- 



2 7 + ,5' ' (7 + 5) 2 ' ' (7 + 5) 

B,A _ A 1 + 7^ , \ 9,A_ l + 7<> \ 

1/4 - 65(5 + 7 ) +A ' * " 6 5(5 + 7 ) ' 
e,A = (1 + 75)^-7) 1-0 _ (l- 7 ) A 3 5 2 -7 
1 ~ 3(5(7 + 5)2 1 4 J 6(7 + 5) 2 7 + 5' 

e = (l+ 7 5)(5 2 - 7 ) n l-_0, _ (1-7) 

2 ~ 35(7 + 5) 2 [ 4 ' 12(7 + 5)' 

In this paper, we also consider models with (formally) lower order accuracy, as defined below. 

Definition 2.4 (lower order decoupled approximations). Let £°,v° be given scalar functions, and 
set parameters (e, /i,7, 5) € "P, (A,#) G K 2 . ^4 decoupled approximate solution of the system (1.1) is 

t/ = ^+(t,a;-t) + w_(t,a; + t),(7 + 5)(u + (t,a;-t)-v_(t,a; + t))^, 

where v± \t=a = ^(C° + ^pj) ari< ^ f± = (1 i ^Ad 2 ) _1 i>± wit/i v± satisfying a scalar evolution 
equation. In the following, we consider 



• the inviscid Burgers' equation: 

(2.13) d t v x ± e ai vid x vi = 0; 

• the Korteweg-de Vries (or more precisely Benjamin-Bona-Mahony) equation: 

(2.14) d t vi ± eait^Q,.^ ± fiv^d^ - ^v e t ' x d 2 x d t vi = 0; 

• the extended Korteweg-de Vries equation: 

(2.15) dtv^ ± ea lV ^d x v^ ± e 2 a 2 (v±) 2 d x v± ± vu£ x d%v$. - ^v* d 2 x d t v\ = 0; 



where the parameters satisfy the identities of Definition 2.3 



The existence, uniqueness and rigorous justification of such family of approximate solutions of 
system (1.1), is investigated in the following section. 

2.2 Rigorous statements 



In this section, we rigorously justify the decoupled approximations defined in Definition 2.3 and 2.4 



in the sense of consistency. We first prove that the evolution equations we consider are well-posed 
for large times, as they satisfy energy estimates in Sobolev norms. Persistence property of estimates 
in weighted Sobolev norms are also investigated, as spatial decay at infinity {i.e. localization in 
space) plays an important role in the accuracy of the approximations, as we shall see. 

These estimates are valid for long times, and we demand the approximation to be valid for 
such times. In order to deal with this, we introduce (following the decomposition introduced 
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in the previous section) an explicit first order corrector, which roughly reproduces the coupling 
effects between the two counterpropagating waves. Including this correction yields an approximate 
solution, which is precise (in the sense of consistency) at high order. The key ingredient is then 
to obtain precise estimates on this corrector, and especially control its the secular growth (which 
should be sublinear in time). 



Let us first study the well-posedness of the Constantin-Lannes equations (2.12), that we repro- 
duce below: 

(2.16) (1 — fj,pd 2 )d t u + ect\ud x u + e 2 a%u d x u + e 3 a^u d x u 

+ n vd x u + fie d x (nudlu + b(d x u) 2 ) = 0. 



Our result, concerning the existence and uniqueness of strong solutions of (2.161, requires that 
> 0, as it relies heavily on a priori estimates of the solution in the following scaled Sobolev 
norm: 

I 1 2 I 1 2 , , 2 

\u\ „ s+ i = \u\ H3 + /i/?|u|„ 3+1 , for some s > 0. 
We will also use weighted Sobolev norm, defined through the following norm: 



= J2\ xju \m 



\x. 



+2(n-j). 



Proposition 2.5 (Well-posedness). Let u° G H s+1 , with s > sq > 3/2. Let the parameters be such 
that P, fj,,e > 0, and define M > such that 

£+~+/i + e+|ai| + |a 2 | + |a!3| + M + N + N < M. 

Then there exists T = C( sg } 3/2 , M, \ u° | ffS+1 ) and a unique u € C° ( [0, T/e) ; H s + l )nC 1 ([0, T/e) ; H £) 

such that u satisfies ( |2.16| and initial condition u\t—o — M °- 
Moreover, u satisfies the energy estimate for < t < T/e: 

(2.17) \\dtu\\ Lxi[0>T/e) . Hp + ||«|L. ([0ir/e) . H .+i ) < C, (^3372' M 'l u °lff;+ 1 ) 



Assume additionally that for fixed n,k £ N, the function x^u° € H s+S , with < j < n and 
s = k + 1 +2(n — j). Then there exists T = C( Sg 2 3 / 2 i M, n, k, Y^j=o u ° \ H =+ k + 1 + 2 (^-i)) such that 
for < t < T x min(l/e, 1/fJt), one has 

1 n 

\\ xn dkd t u\\ L ^ m) . H ^ + \\x n d k u\\ Lao(m . H3 + t) < g( — 3/2 , M, n, k, g \x j u°\ K+k+1+2 ^-^ . 
Ln particular, one has, for < t < T x min(l/e, 1/fX), 

(2.18) ||Hl£- ([ o, t)i x SiM ) + ll«IL«-( [ o,t) s ^) ^ C (^372' M '"'I W °I^)- 

Remark 2.6. A^oie i/iat </ie above result allows the parameters (except jifi) to vanish, so that 
we ensure the well-posedness and control in Sobolev norms of all our decoupled models defined in 
Definition\2.3\ and\24\ 

Of course, such results are already well-known for inviscid Burgers', KdV and eKdV equations 
(see for example ]1(A \3R [23 \2% and references therein), and actually do not require fi/3 > 0. 
The case of Constantin-Lannes equation is given in Proposition 4 in JJBj- 

The persistence of the solution in weighted Sobolev norms for the Camassa-Lannes equation is 
new, as far as we know. Similar results in the case of (eventually extended) Korteweg-de Vries equa- 
tions are obtained in slightly different setting (for the most part using weighted L 2 spaces intersected 
with non-weighted Sobolev spaces H s , s > 0) in JM\M\M\M\U\\Z^\M U3- 
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Propositi on|2.5| is proved in Appcndix[B] We now state our main result, whose proof is postponed 
to Appendix [C| 

Proposition 2.7 (Consistency). Let (°,v° £ H s+e , with s > s > 3/2. For (e,ju,5,7) = p £ V, 
as defined in (1.2), we denote Uq l the unique solution of the CL approximation, as defined in 
Definition 2.3 For some given M* +e > 0, sufficiently big, we assume that there exists T* > and 
a family (Uqi^p&T such that 

T* = max(T>0 such that \\U^ L \\ LOS{[QT) . HS+<i) + \\dtU% L \\ Ltx>{[0>T) . HS+5) < M* +6 ). 

Then there exists U c — U C [U^, L ] such that U = Uq L + U c is consistent with Green-Naghdi 
equations (1.1) of order s on [0,t] for t < T* , at precision 0{e* CL ) with 

e* CL = Cmax( f 2 (^ 7 ) 2 , f V 2 ) (l + Vi), 
with C — C( sq _ 3 , 2 , M* +6 , j-'— , <5 max , e max , /i max , |A|, \8\), and the corrector term U c is estimated as 

II^IL-dCT*];^) + ll 9tC/C |L~([o,T*];_f/o - C max(e((5 2 -7),e 2 ,^)min(t,v / t). 



Additionally, if there exists a > 1/2, AfjLg, T" > such that 

6 5 

£ Ha +^ 2 r^^ i || ioo([0iT);ff3) + ^ ||(i +x 2 ) a d%u» CL \\ LOO([QT) . HS) < m« +6 



k=0 k=0 

then U = Uq l + U c is consistent with Green-Naghdi equations (1.1) of order s on [0,t] for t < T*, 
at precision 0{e\j L ) with 

4 £ = C max(6 2 ( ( 5 2 - 7 ) 2 ,6 4 , A1 2 ), 
with C = C( SQ _ 1 3/2 ,Mf +6 , fi^,5 max , e max , /Lt max , |A|, \0\) and U c is uniformly estimated as 

\\ UC \\l°°(10,t»];H>) + ll 5 * C/C |lz,~([o,rtt];ij-) - C max(e((5 2 - 7), e 2 , min(t, 1). 

Remark 2.8. The function U c — U c [Ucl] is given explicitly in the proof of the Proposition; see 
Appendix\C\ and more precisely Definition \Gl\ 

Remark 2.9. In the estimates presented in Proposition \2. 7| and Cor pilar y \2.1C\ below, the terms 
growing in 0(yt) come from coupling effects between the two propagating waves, that are neglected in 
our decoupled models; this is why the accuracy is significantly better if the initial data is sufficiently 
spatially localized, as the two counterpropagating waves will be located far away from each other 
after long times. 

The function U c , which depends only on the decoupled waves v±, is a first order corrector 
which allows to take into account the leading order coupling effects, and therefore reach the desired 
accuracy. 

Uniformly bounded terms are the contribution of unidirectional errors, generated by the differ- 
ent manipulations on the equation (e.g. BBM trick), and eventual neglected terms in lower order 
approximations in Corollary \2.1C\ The magnitude of each contribution depends on the situation 
(small amplitude of deformation, shallow water, critical ratio, etc.) and is discussed for different 
scenari in subsection\2.3\ below. 



Equivalent results may be obtained for lower order approximations, when following the steps of 
the strategy presented in Appendix [B] More precisely, this modifies mostly the last term of ( C.5 ) in 



Lemma |C.3| adding the contribution of terms that are neglected in the chosen evolution equation. 
The additional error is therefore uniformly bounded on times [0, T*] and [0,T*]. The Corollary 
below details the accuracy of such approximations, and its proof is omitted. 
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Corollary 2.10. Assume that the hypotheses of Proposition 2.7 hold. Denote U^ Kdv , U v Kdv and 



U[ B , respectively, the solutions of the eKdV, KdV and iB approximations, as defined in Defini- 



tion 2.4 In each case, we assume that the decoupled approximation is uniformly estimated in 



[0, T*], as in Proposition \2. 7| Then 
i. there exists U c = U c [U^ Kdv ] such that U = U^ Kdv + U c is consistent with Green-Naghdi 



equations ( 1.1 1 of order s on [0, t] for t < T* , at precision 0(e* Kdv ), with 

= C x ( max(e 2 (<S 2 -7) 2 ,e 4 ,Ai 2 ) (1 + Vt) + max(e 3 ,^e)) 



-eKdV 



ii. there exists U c = U c [Uxdv^ suc h that U = U v Kdv + U c is consistent with Green-Naghdi equa- 
tions (1.1) of order s on [0,t] for t < T* , at precision 0{e* Kdv ) with 



£ KdV 



C x ( max(e 2 (£ 2 - 7 ) 2 ,e 4 , M 2 ) (1 + yft) + e 2 ) ; 



Hi. there exists U c = U c [Uf B ] such that U = Uf B + U c is consistent with Green-Naghdi equa- 
tions (1.1) of order s on [0,t] for t < T* , at precision 0{e* B ) with 



C X ( max(e 2 ( ( 5 2 - 7 ) 2 ,e 4 ,^ 2 ) (1 + yft) + max(e 2 ,^)) ; 



where C = C(^^, M; +6 , j^—, S max , e max , /i max , |A|, \6\). Each time, the corrector term U c is 
estimated as follows: 

ll ?/C |L~([ci,T*];ff') + ll 5tC/ 1L-([o : T*];H=) - C max(e(<5 2 -7),e 2 ,/i)min(t,^). 
Moreover, if the uncoupled approximation is suffic ientl y localized in space, then all the estimates 



are improved as in the second part of the Proposition 2.1 (that is replacing yft by 1). 



Note that Proposition 2.5 ensures that Proposition |2.7| and Corollary |2.10| are not empty, but 



on the contrary are valid for long times, provided that v t ' X > and the initial data sufficiently 



smooth. More precisely, the following result is a straightforward consequence of Proposition |2.5| 
Corollary 2.11 (Existence and magnitude of T* , T«). Let ((°,v°) = U° £ H s+7 ,s > s > 3/2, 

A 

and let p £ V , with additional restriction v t ' > vq > 0. 

Then there exists C ll C 2 = C( S(t 2 3 / 2 > /'max, e max , l/S min ,S max , |A|, \9\,l/u , \ U°\ H s^+7), indepen- 
dent of p, such that for any M* +6 > C\, the decoupled approximate solutions defined in Proposi- 
tion 2. ? and Corollary \2.10 are uniquely defined and satisfy the uniform bound of the Proposition, 



with 

T* > C 2 /e. 

If U° £ X?, +1 , then there exists Ci,C 2 = C( So ^ 3/2 , fam^mmCn^max, |A], l^l,^ 1 , |[/°| x , +7 ) 
such that for any M* +6 > G\ , one has 

T« > C 2 / max(e,/i). 



One question, which is essential in the discussion of section |2.3[ is whether these above estimates 
are optimal. Let us discuss some elements of answer. First, it is well-known that inviscid Burgers' 
equation d t u + ud x u will generate a shock in finite time, for any non-trivial, decreasing at infinity 



initial data. A simple scaling arguments shows that (2.13) will therefore generates shock in finite 
time T s=! 1/e. Adding a linear dispersive term as in the KdV equation regularizes the solution, 
and in particular induces global well-posedness for sufficiently smooth initial data. Moreover, there 
exists an infinite number of conserved quantities [471 148] , which allow to control uniformly in time 
the Sobolev norm of the solution ||u|| ffs , s £ N. This result is true for the whole class of extended 

Korteweg-de Vries equations (2.15), if vp X = 0. 



On the contrary, it is known that Camassa-Holm family of equations, related to (2.12) can 



develop singularities in finite time in the form of wave breaking |18] . In |19j , it is shown that wave 
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breaking occurs for the solution of (2.12 1 for a specific set of parameters, in time T ss 1/e, provided 
that the initial data is sufficiently big in L°° norm. However, as the justification of our model 
assumes that the initial deformation is bounded in Sobolev norm, uniformly with respect to the 
parameters, assumptions of [19j Proposition 6] cannot be justified. We also let the reader refer 
to [12] , and references therein, for an in-depth study of the conditions (in term of decay at infinity 
of the initial data) for finite time blowup of the solution of Camassa-Holm equation. 



2.3 Discussion 



Our result takes its full meaning in the light of a conjectured stability result on Green-Naghdi 



equations (1.1), or any consistent model, that would have the following formulation: 



Hypothesis 2.12 (Stability). For (e, /i, S, 7) = p € V , as defined in (1.2), let Let Uq N be a solution 



of the Green-Naghdi system (1.1) such that (Uq N ) is uniformly bounded on H s , s sufficiently big, 



over time interval [0,T]. Then let (U v ), satisfy U v \t—o — UQ N \t=o, be consistent with Green- 
Naghdi equations (1.1) on [0, T], of order s at precision 0(e) (see definitional^. Then the difference 
between the two family of functions is estimated as 



with C = C( sup„ 



\U P I 

\ u gn\ 



\U p 



L°°([0,T];H') 



U P I 

U GN\ 



L°°([0,t];H° 



: C e t, 



sup„ \\U p \ 



L°°([0,T];ff=+ 5 )' <5„ 



/in 



Such a result has been obtained for a well-chosen Green-Naghdi system in the case of the water- 
wave problem in |441 l2l 134] . and is under current investigation by the author and collaborators. 

If the stability Hypothesis |2.12 holds, one can deduces from the consistency result of Proposi- 
tion 2.7 that the difference between the solution of the Green-Naghdi system ( 1.1 ) (and in the same 



way, the solution of the full Euler system if it exists; see Remark A. 6) and the (weakly coupled) 
approximate solution U = Uql + U c is small. The estimates concerning the difference between the 
solution of ( |1.1[) a nd the fully decoupled solution Uql simply follows from the estimates on U c , 



in Proposition 2.7 This strategy has been used in the water-wave case in [7], and in the case of 



internal waves (when restricting to the long wave regime) in |23j . 

Throughout this section, we assume that Hypothesis |2.12| holds, and study the convergence 
results between solutions of the Green-Naghdi model and the different approximate solutions which 
proceeds. 



Let us first state the convergence results concerning coupled models. 



Corollary 2.13 (Convergence between coupled models). For (e, /i, 5, 7) = p € V, as defined 
in (1.2), let Uq N be a solution of Green-Naghdi equations (1.1) such that the family (Uq N ) is uni- 



formly bounded on H s , s sufficiently big, over time interval [0,Tgn]- Let f/| = (u p +u p , (-f+5)(u p — 
uf,)) whereu p ,u p are the solutions of coupled equations (2.3a)-(2.3b) such thatU c \t—o = Uqn \t=o , 
and assume that the hypotheses of Proposition 2.2 hold over time interval [0, T c ]. Then if Hypothe- 
sis 



2.12 is valid, one has for any t < mm(TQN, T c ): 



\^GN 



L°°([0,t];H°) 



< C (e 4 + fie 2 )t , 



with C = C(sup p ||C/£ 



GN\\L°°([0,T GN ];H°)> SUP P Wc I L°°([0,T c ] <5 mi „ '^max, tmax^max) 



Proof. This is a direct application of the consistency Proposition 2.2 together with the stability 
Hypothesis [2~T2l □ 



Remark 2.14. Note that the consistency result of Proposition pO| allows to obtain a similar es- 
timate on the difference between the solution of the Green-Naghdi system (1.1) and the solution 
of the full Euler system (A.4) with same initial data, if the latter exists (see the 
Remark A. 6). This estimate is of size 0(fi 2 t), 



if the latter exists (see the discussion in 
so that triangular inequalities allow to extend the 
estimates of Corollary \2.15\ and \2.lb\ below, replacing Us by the solution of the Full Euler system. 



Let us now turn to the weakly coupled model, defined in Proposition 



2.7 
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Corollary 2.15 (Convergence of weakly coupled model). For (e,fi,5,j) = p e V, as defined 
in (1.2 1, let U GN be a solution of Green-Naghdi equations (1.1) such that the family (U GN ) is 
uniformly bou nded on H s , s sufficiently big, over time interval [0,Tqn]. Assume that hypotheses 



of Proposition 2.7\ hold, and denote U v = U GJ + U C [U GL ] the approximate solution, satisfying 



Ugn p |t=o = U v [t=o ■ Then if Hypothesis 2.12 is valid, for any t < mm(T G N,T* +6 ), one has 



with C = C(su Pp \\U GN \\ Loo(l0TGN] . Hsy Mt 



|A|,|0|). 



P II^GW llL°°([o,T G N];-fi" s )' """ s+6 ' J ™ ' max ' £max > 

If moreover, the initial data is sufficiently localized in space, then one has for t < mm(T G N, T^ +6 ), 



\ttP_tjP II 



< C max(e 2 ((5 2 - 7 ) 2 ,e 4 ,^ 2 ) t, 



s-\-Qi 5 m in ' Umax ' c niaxi /^max 



c|A|,|0|). 



with C = C(su Pp \\U GN \\ LBca0!T] . H3y M : 

We therefore see that the weakly coupled model may achieve the same accuracy as the fully 
coupled model, in the critical case, and if the initial data is sufficiently localized in space. 

Let us now turn to the fully decoupled models. The following result is a straightforward appli- 
cation of the above Corollary, together with the estimate of Proposition 2.7 Estimates concerning 



lower order decoupled models are obtained in the same way, using Proposition 2.10 



Corollary 2.16 (Convergence of uncoupled models). For (e, /x,£,7) = p G V . as defined in (1.2), let 
U GN be a solution of Green-Naghdi equations (1.1) such that the family (U GN ) is uniformly bounded 
on H s , s sufficiently big, over time interval [0,Tgn]- Denote respectively U GL ,U^ Kdv ,Ux dv ,U[ B 
the decoupled approximations defined in Definitions \2.3\ and \2.4\ Assume that the hypotheses of 
Proposition 2.7 hold, and Hypothesis 2.12 is valid. Define So — max(e((5 2 — 7),e 2 ,/i). Then for any 
t < min(T(3iv,T* +6 ), one has 

L~([0,t];i? s 



\u CL -u GN *\ 



\u v 

\ u eKdV 

\\U V 
\\ u KdV 

\\U P 
\\ u iB 



U GN f 
U GN f 
U GN f 



lL°°([0,t];ff s 
L=°([0,t];ff s 



lL°°([0,t];ff s ) 



< C s min(t,i 1/2 )(l + e t), 

< C £ min(t,i 1/2 )(l + e t) + C max(e 3 , fxe) t, 

< C e min(t,i 1/2 )(l + e t) + C e 2 t, 

< C s min(M 1/2 )(l + s t) + C max(e , fi) t, 



with C — C{ sup p \\U GN \\ L ^^ TcN y Hs yM* +6 , g^— , <5 max , e max , ^i max , |A|, \9\) . 



If the initial data is sufficiently localized in space, then one has for t < min(T, T"), 
\\ u CL-U s \\ LOO([o t] . HS) < C so min(t,l)(l + e t), 

\\UeKdV 



U s\\ L °°{lo,t);H>) < C £ o min(t,l)(l 



\Ukav — Us\ 



iB 



Us 



£°°([0,i];H*) 
i«([0,t];H') 



with C — C[ sup 



P l\ U GN\\L™([0,T GN ]■,H 3 ) , ' 



< C so min(t, 1)(1 

< C s min(t, 1)(1 



s t) 
s t) 
s t) 



C max(e , fie) t, 



C 



t. 



s _|_6? 5 min ' u max: ^maxj Pmax? 



I- C max(e 2 ,//) t, 
|A|,|0|). 



Corollary |2. 16| exhibits different sources of error with different time scales. 

• The first (and often main) source of error comes from the coupling between the two counter- 
propagating waves which are neglected by our decoupled models, but recovered at first order 
by the weakly coupled model (see Corollary 2.15). Following [52] . we name the contribution 
of this term the counterpropagation error. 

— This error grows linearly in time, for times of order 0(1), as the two waves are located 
at the same position; and coupling effects are strong. 

— However, one is able to control sublinearly the secular growth of such coupling terms, 
which is the key ingredient of our result, and yields a contribution of size O(soVi) in general, 
O(so) if the initial data is sufficiently localized. 
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residual 
coupling 
effects 




(a) non-localized initial data (b) localized initial data 

Figure 1: Sketch of the error 



• At very long time, one sees the effect of the precision of the consistency result, through the sta- 
bility hypothesis. Such a contribution is unavoidable, as it appears after many manipulations 
of the equations, such as the use of BBM trick, or near-identity change of variables. We call 
these contributions residual errors. The error generated in that way affects both the solution 
of the scalar evolution equation and the coupling term, which creates two contributions. 

- The so-called unidirectional error is linear in time as long as the solution of the 
decoupled model is uniformly bounded, and varies greatly with the choice of the model. 

— The residual error coming from the coupling term may be superlinear (if the initial 
data is not localized in space), as the coupling term grows in time. 

All of these contributions are summarized in Figure [I] The total error of a decoupled models is the 
sum of the green and the red curves. The green curve represents the error generated by neglecting 
the coupling between two waves, and therefore is affected by three of the contributions listed above, 
at different time scales. The red curve is the unidirectional contribution, and varies with the choice 
of the model. 

However, let us precise that the full pattern is not always visible, as in many cases, one source 
of error will be negligible in front of the others. In particular, keep in mind that we obtained 
existence and uniform control in weighted Sobolev spaces of the solutions of decoupled models only 



over times of order 0(1/ max(e, /i)) (Proposition 2.5|, thus the time where residual contributions 
appear is always out of the scope of our rigorous results. 

In order to understand the relevance of a specific evolution equation, one has to look whether 
the unidirectional contribution, which depends on the evolution equation itself, is smaller or greater 
than the different sources of error due to coupling. For example, unidirectional error is always 
negligible in front of coupling terms in the Constantin-Lannes model, whereas it is expected to be 
preponderant in the inviscid Burgers' model. 

Let us look precisely at several interesting scenari, in the following subsections. Each discussion 
is supported by numerical simulations. Each time, we compute the coupled Green-Naghdi model 
and various evolution equations — among iB, KdV, eKdV, CL — , for different values of e and 
corresponding /j, (// = e or fj, = ^Je, depending on the situation), and over times C(e _3/2 ). As 
pointed out above, such time scale is out of the scope of our rigorous result. However, interesting 
behavior appears after times 0(l/e), as we shall see. Each of the figures below contains three 
panels. On the left-hand side, we represent the difference between the Green-Naghdi model and 
decoupled models, with respect to time and for e = 0.1,0.05,0.035. Values at times 1/e and 1/e 3 / 2 
are marked. In the two right-hand-side panels, we plot the difference in a log-log scale for several 
values of e (the markers reveal the positions which have been simulated), at given times 1/e and 
1/e 3 / 2 . The pink triangles express the convergence rate. 

The initial data is fixed such that the left-going wave is initially two-third the right going wave, 
to avoid symmetry cancellations: i>_ | t=0 = 2/3i> + \t=o- For localized initial data, we choose initial 
data v + \ t= o — exp(— (x/2) 2 ) and for non-localized initial data, we choose v + \ t= o = (1 + 10a; 2 ) -1 / 3 . 
We set 5 2 = 7 = 0.64 in the critical ratio setting, and 8 = 0.5,7 — 0-9 i n the non-critical ratio 
setting. Finally, we set 9 = 1/2 and A = 0. 
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The major difficulty of the numerical simulations computed throughout this paper is the fact 
that as e is small (down to 0.035 in our simulations), the time domain t G [0, 1/e 3 / 2 ], and therefore 
the space domain (as the two waves are moving at velocity c± ~ ±1) of computation becomes very 
large. The decoupled models can be solved very efficiently by using a frame of reference moving 
with the decoupled wave, but solving numerically the Green-Naghdi scheme has to be carried out on 
the full time/space domain. Thus we need a scheme which allows a great accuracy, with a relatively 
small computation cost. With this in mind, we turn to multi-step, explicit and spectral methods. 
The space discretization, and in particular the discrete differentiation matrices use trigonometric 
polynomial on an equispaced grid, as described in [59] (precisely (1.5)). This yields an exponential 
accuracy with the size of the grid Ax, if the signal is smooth (note that the major drawback is 
that the discrete differentiation matrices are not sparse). It turns out setting Ax = 0.2 is sufficient 
for the numerical errors to be undetectable in our computations. After this space discretization, 
one has to solve a system of ordinary differential equations in time, and we use the Matlab solver 
odell3, which is based on the explicit, multistep, Adams-Bashforth-Moulton method [57], with a 
stringent tolerance of 10 -8 . 

In the following discussion, we denote by Co a constant, which do not depend on all the other 
parameters involved, and we write A = 0(B) if A < C Q B, and A w B if A = O(B) and B = 0(A). 



2.3.1 The long wave regime. 



Let us assume that we are in the long wave regime: e = O(fi), and therefore £q ~ \i in Corollary 2.16 



It follows that the decoupled KdV approximation has the same order of accuracy as the higher 
order models, whatever the initial data or critical ratio is. Indeed, one has the same estimates for 
the decoupled approximations U = Uxdv or U — U e Kdv or U = Ucl- 

\\ U - U GN\\ L o O([0)t] , H , ) < Con min(i,i 1 / 2 )(l + id), 

for any t € [0,T* +6 ), and if the initial data is localized, 

\\ u - u GN\\ Laoi[o t] . Hs) < C n min(*,l)(l + /it), 

for any t £ [0, T s " +6 ). Let us note that one recovers the results of [53], for the KdV approximation. 



E = 0.1 0.05 0.035 








0.035 0.05 0.075 0.1 
epsilon 



0.035 0.05 0.075 0.1 
epsilon 



(a) behavior of the error with respect to time 



(b) behavior of the error with respect to e = fi 



Figure 2: Long wave regime, critical ratio, localized initial data 



More precisely, a unidirectional error of size 0(ii 2 t) is not seen, as it is smaller than the coupling 
terms, presented above, so that the additional error produced by neglecting higher order terms in the 
Constantin-Lannes approximation [2~12] and considering only the KdV equation, does not change 
the accuracy of the model. However, let us note that for localized initial data, unidirectional error of 
size 0(jjL 2 t) is the limiting case. This is why one sees in Figure 2(a) a noticeable difference between 
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the K dV ap proximation and higher order approximations (eKdV, CL). The rate of convergence in 
figure 2(b) however, is identical: the error is of size O(n) at time T = 0(1/ n) and 0(n 1/2 ) at time 
T = 0(l//.t 3 / 2 ). We remark that the last panel of Figure [2] shows a slight discrepancy with respect 
to the predicted convergence rate in 0(fi 1 ^ 2 )- This may be due to the fact that higher order sources 
of error (typically of size 0(fj. 3 t)) are detectable for larger values of fi (/x = 0.1), and eventually 
becomes negligible for smaller values (fi = 0.035), therefore artificially improving the convergence 
rate. Let us note also that the criticality of the depth ratio do not play a role in this analysis, and 
that simulations in the case of non-critical ratio gives similar outcome. 



2.16 



2.3.2 The Camassa-Holm regime, with non-critical ratio. 

We are now in the case where \S 2 — j\ > a Q > 0, and e ^//x, thus £0 ~ e ~ v/M m Corollary 

In that case, the contribution of the coupling error are always greater than the one of the 
unidirectional error, and one has the same estimates as for the decoupled approximations U = UiB 

U = U K dV, U = U e KdV or U = U C L- 



U-U GN 



L°°i[0,t);H*) 



< C y/ji min(t,t 1/2 )(l + y/jit), 



for any t £ [0,T*, 6 ), and if the initial data is localized, 



\U-Ugn\ 



L°°([0,t];Hn 



< C yf]i min(t,l)(l + y/jit), 



for any t £ [0,Tg, 6 ). As a consequence, the inviscid Burgers' approximation is as precise as any 
higher order decoupled model. 



e = 0.1 0.05 0.035 





I 10"' 



0.035 0.05 0.075 0.1 
epsilon 



0.035 0.05 0.075 0.1 
epsilon 



(a) behavior of the error with respect to time (b) behavior of the error with respect to e = ^fjl 

Figure 3: Camassa-Holm regime, non-critical ratio, localized initial data 

Here, the unidirectional error of the inviscid Burgers' approximation is of the same order of 
magnitude as the residual coupling error, if the initial data is sufficiently localized. Therefore, in 
figure [3j one can observe a noticeable difference between the iB approximation and higher order 
models for long times, although the convergence rate is similar at time T = 1/e = 1/y/ji. For longer 
times, the decoupled approximate solutions do not seem to converge, which may indicate that the 
exact solution Uqn cannot be controlled over times 0(e~ 3 / 2 ), in the setting we use. 



2.3.3 The Camassa-Holm regime, critical case. 

Now let us assume that 5 2 — 7 = O(JJi), and e w JJi, so that £0 



jji in Corollary 



2.16 



We recall that T* and T" are known to exist and to be of size T = 0(1/ max(e, /j,)) = 0((i 1 / 2 ), 



as a result of Proposition 2.5 Over such time, the unidirectional error in the eKdV model is smaller 
than the coupling errors: one has for both U = U e Kdv and U = Ucl, 



U-U GN 



Z,°°([0,t];ff») 



< C H min(t,i 1/2 ), 
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for any t £ [0,T* +6 ) if T* — 0(fi 1 / 2 ). If moreover, the initial data is localized, then 

||^ -C7 Gw|L=o([o,t] ; ff») ^ C n min(t,l), 

for any t € [0, r| +6 ) if T" = C(/i _1 / 2 ). The accuracy of both eKdV and CL approximations are 
p = 0(> 3/4 ) at time T = 0{^ 1/2 ) = C(e _1 ), or p = 0(n) if the initial data is sufficiently localized. 
The accuracy of the KdV and iB approximation is ©(/i 1 / 2 ) at that time. 

However, if one looks at longer times, then the picture is different. In our simulations, we looked 
at times up to T = 0(e~ 3 / 2 ) = 0(fi^ 3/>4 ). At that time, the localization in space of the initial 
data plays an important role. If the initial data is non-localized in space, then contribution of the 
coupling terms are predicted to be greater than the unidirectional error of the eKdV approximation: 
the secular coupling error is dominant up to T = 0(fi^ 1 ) (/J.\/t > /i 3 / 2 i), after what the residual 
coupling term is dominant (/_t 2 i 3 / 2 > p?l 2 t for t > j-i^ 1 )- Thus the eKdV approximation is predicted 
to be as precise as the CL approximation. On the contrary, if the initial data is localized, then 
the unidirectional error is dominant: /i 3 / 2 t > + fit) for any t > ^T 1 / 2 . This leads, at time 
T = 0(/i~ 3 / 4 ), to an error of size 0(/j,) for the CL approximation, and of size C(/j, 3 / 4 ) for the eKdV 
approximation. This phenomenon is clearly seen when comparing Figure [4] (localized initial data) 
and Figure (non- localized initial data). 




50 100 150 0.035 0.05 0.075 0.1 0.035 0.05 0.075 0.1 

time epsilon epsilon 



(a) behavior of the error with respect to time (b) behavior of the error with respect to e = ^/JL 

Figure 4: Camassa-Holm regime, critical ratio, localized initial data 



E = 0.1 0.05 0.035 T=0fe T=0/e : 




50 100 150 0.035 0.05 0.075 0.1 0.035 0.05 0.075 0.1 

time epsilon epsilon 



(a) behavior of the error with respect to time (b) behavior of the error with respect to e = ^/JL 

Figure 5: Camassa-Holm regime, critical ratio, non-localized initial data 
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3 Unidirectional propagation 

In the previous section, we proved that approximate solutions of the internal wave problem consists 
in two decoupled, counterpropagating waves, each of them evolving according to a given scalar 
evolution equation (Definitions 2.3 and 2.4). However, it has been seen that the accuracy of such 
approximation may be low, especially in the Camassa-Holm regime e ~ yfji, if the ratio is non- 
critical (5 2 ^ 7) and/or the initial perturbation is not spatially localized (see section 2.3.2), as 
significant coupling effects develop between the two counterpropagating waves. 

In this section, we show that if one chooses carefully the initial perturbation (deformation of 
the interface as well as shear layer-mean velocity) in order to root out one of the two waves , then 
the accuracy of the uncoupled model (that we call unidirectional approximation in that case) can 
be considerably improved. This study follows the strategy developed for the water-wave problem 
in PHE5]. 

The result is precisely stated below, followed by its proof and a brief discussion. 

Proposition 3.1. Set A, 9 £ R, and (° £ H s+5 with s > s > 3/2. For (e,//,<5, 7 ) = p £ V, as 
defined in (1.2), denote (C v )pev ^ e unique solution of the equation 

(3.1) d t ( + d x C + eaxCdxC + e 2 a 2 C 2 d x ( + e 3 a 3 ( 3 d x C + AU# A flgC - /wf'^ftC 



with 



3 5 2 - 7 _ 21(5 2 - 7 ) 2 5 3 + 7 

ai ~ 2^+7' " 2 ~ 8(7 + 5) 2 7 + <5' 

_ 71(5 2 - 7 ) 3 37(5 2 - 7 )(5 3 + 7 ) 5(5 4 - 7) 

16(7 + 5) 3 4(7 + 5) 2 + 7 + 5 ' 

^Hl-^A)^, ^ A = (« + A)- 1 + 7 * 



6(5(7 + 5)' t v ■ ^5(7 + 5)' 

j,x = (14-6(fl + A))(5 2 - 7 )(l+ 7 <5) 1-7 

n 245(7 + S) 2 6( 7 + 5) ' 

6iX _ (17-i2#)(5 2 - 7 )(l + 7 <5) I-7 



2 48<5( 7 + c5) 2 12(7 + 5)' 

For given M s+5 , h > 0, assume that there exists T s+5 > smc/i i/ia< 

T s+5 = max(T>0 such that ||C P || L oo( [0>T ). H a+5) < ) , 

and for any (t,x) £ [0,T s+ s) x R, 

fti(t,af) = 1 - eC P (^^) >h>0, h 2 (t,x) = \ + eC P (^^) > to > 0. 

5 

TTien de/me v? as v<> = hl +£ 2 v[C], with 

(3.2) v[Q = C + £ yC 2 + e 2 yC 3 + £ 3 f C 4 + M^ 2 C + M £ (^C + K2 (^C) 2 ), 

where parameters 0*1,0:2,0:3 are as above, and 

l + 7 <5 14(5 2 - 7 )(1 + 7 5) 1-7 17(5 2 - 7 )(l+7<5) I-7 

"~ K l = ?T777 1 cVo 7T, — ni K 2 ~ 



65(7 + 5)' 245(7 + 5? 6(7 + 5) ' 485( 7 + 5) 2 12( 7 + 5) ' 



Then {C^tV^) is consistent with Green-Naghdi equations (1.1), of order s and on [0,T s+ 5) ; with 
precision 0(e), wii/i 

£ = C(M s+5 , ^— — — *— , 5 max ,e max ,/i max , |A|, |6»|) x max(e 4 ,/i 2 ). 

S " 3/2 5 mirl 

Remark 3.2. When 7 — > and 5 — > 1, one recovers the one-fluid model presented in JTPI sec- 
tion #.,2/, luii/i a = -^g^ and A = 0, using notations therein. 



20 Decoupled and unidirectional asymptotic models for the propagation of internal waves 



Remark 3.3. In Proposition \3 . 1\ the approximation consists in a scalar evolution equation satisfied 
by the deformation at the interface, and we reconstruct the shear layer-mean velocity from the defor- 
mation (and in particular, this determines the initial shear velocity, from the initial deformation). 
Following [10], a similar strategy could consist in looking for an evolution equation for the shear 
layer-mean velocity, and reconstruct the deformation at the interface. We decide not to present the 
outcome of such strategy, as the result is very similar, and calculations somewhat heavier in that 
case. 



Remark 3.4. As discussed in \ 191 Proposition 5], specific values of parameters in (3.1| yield 
equations with different properties, especially concerning the behavior near the maximal time of 
definition (if it is finite). Indeed, the proof of f!9l Proposition 5] can easily be adapted to more 
general coefficients, and one obtains: 

• If v t ' > 0, k x ' = 2k2 ,A > and > 0, then singularities can develop in finite time only in 
the form of surging wave breaking. In other words, if the maximal time of existence of £ is 
finite, T < oo, then 

sup {\C(t, x)\} < oo and sup{d x ((t, x)} t oo ast^T. 

te[Q,T),x€R i£l 

• If Vt' X > 0, Ki' X — 2k2 ,A < and a 3 < 0, then singularities can develop in finite time only in 
the form of plunging wave breaking. In other words, if the maximal time of existence of £ is 
finite, T < oo, then 

sup {|£(i, x)|} < oo and inf {d x ((t, x)} J, — oo ast^T. 

te[o,T),xeR zeR 

Identity K^' — holds in the line 9 — A = 1/2, and in that case, u t ' > if and only if 

9 > 1/4. Restricting to 9 < 1 as natural values for the use of BBM trick, one can easily check 
that if 7 = 0, then singularities may occur only as surging wave breaking, as it is the case in the 
one-layer situation. On the contrary, 2/7 ~ 1, as is the case in the stratified ocean (small variation 
of densities) then singularities will occur as surging wave breaking if S > 1 (thicker upper layer), 
and plunging wave breaking will occur for S < 1 (thicker lower layer). 



Proof of Proposition \3.1\ In order to simplify the calculations, we use the Green-Naghdi system ( fTTTj ) 

expressed using the variables (C,V.) where we define v — ^ + \ v. The system reads 

(3.3) 

r dt( + d x v = 0, 



di 



with the following operators: 

Q[hiMv = ^^ 2 {^ A { h ^ d ^)) + 7^ 2 a E (Vc>4^) 

Tl[hiMV = \(( h 2d x (^)) 2 -l(hid x (^) 



Using this system simplifies considerably the analysis. Indeed, it is clear that if £ is to satisfy 
the following scalar evolution equation, 

dtC + d x ( + e ai (d x ( + e 2 a 2 t 2 d x ( + e 3 a 3 ( 3 d x ( 
(3.4) +^d 3 x C + tied x («iC^C + K 2 (d x () 2 ) = 0, 
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then v shall satisfy (using the first equation of (3.3 1, and the fact that the system is at rest at 
infinity: w — >■ when x ~ > ±oo) 

(3.5) v = C + ^^ + ^f C 3 + e 3 f C 4 + ^C + ^(^C + «2(^C) 2 ). 



Now we will show that one can choose coefficients ai,a 2 , etc. such that the second equation of (3.3) 
is satisfied up to a small remainder. 

Indeed, when plugging (|3.4[) and (|3.5|), and expanding in terms of e and \x, one obtains 



- e 2 (6(<5 3 + 7 ) - 5(<5 2 - 7 )a a + (a, 2 + 2a 2 )( 7 + 6))( 2 d x ( 

+ e 3 (l0((5 4 - 7 ) - 9{S 3 + 7 )ai + (2 Ql 2 + 14/3 a 2 )(S 2 - 7) - 2(a 3 + aia 2 )(<5 + l))( 3 d x ( 
,l+7<5 



36 



-2v x (6 + 1 ))d x ( 



+ Lie(4(5 2 - j)u x - 2(a x v + k ± )(S + 7) - 



1-7 2(1+7(5) 



3<5 



(3.6) 



+ ^e(2((5 2 - 7)1/3 - (3ai^ + 2ki + 4k 2 )(<5 + 7) 
= K 



2(1-7) , 2(1+7(5) 



ai)5 x C3 2 C 



where the remainder 1Z can be estimated, provided hi, h 2 > h > 0, as 



\n\ H . < c(\c\ HS+5 , 



1 



s - 3/2 



h~\ 



1 



' 1 "max; ^maxi Atmax) x max(e , fj, ) 



The left hand side of (3.6 1 vanishes when we set 

3 (5 2 - 7 



Oil 



a-2 



5(6 2 - 7 )ai -6((5 3 + 7) _ a| 

2~ 



«3 



27 + r 2( 7 + <5) 

10(<5 4 - 7) - 9((5 3 + 7 )ai + (14/3 a 2 + 2a 2 )(<5 2 - 7) 



2(7 + 5) 
4(<5 2 - 7 H 



1-7 



+ oti 



2(1+7(5) 



35 



K-2 



(S 2 - i)v x 



2(7 + 5) 



U\OL2, V 



U\V, 



1 + j5 
65(7 + S) 



1-7 



1+75 



2(7 + <5) 



2>OLlV Ki 

~ 4 ~2~ 



Such a choice corresponds to parameters of Proposition 3.1 with 9 = A = 0. 



The cases # ^ and A ^ are obtained as in section 2.1 using BBM trick and near- identity 
change of variables. We detail the calculations below, using for simplicity the notation O s (s) for 
any term bounded by £C(|c| ffs+s ). 

BBM trick. We make use of the first order approximation in (|3.4| : 



d t ( + d x Q + eaiC&cC = 03(max(^, e 2 )), 
so that one has, for any 9 £ K, 

d x ( = (l-9)d x (-9(d t ( + eaiCd x () + 3 (max(^, e 2 )). 



Plugging this identity into (3.4 1 yields 
dtC + d x ( + eai(d x C + t 2 a 2 C, 2 d x C + e 3 a 3 C 3 d*C + M(l - 9)v%t ~ ^d 2 x d t C 



(3.7) 



+Hed x ((m - 9v ai )(d 2 x ( + (/t2 - 0i/ai)(0 x C) 2 ) - 5 (max(/i 2 , M e 2 )), 



Conversely, £0, a solution of ( |3.7[ ) (with zero on the right hand side), satisfies (3.4) with a rea- 
minder bounded by 05(max(/i 2 , /ie 2 )). One can easily check that, defining v g as a function of £0 
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3.1 



is 



through (3.5|, {Ce^Ug) satisfies (3.6) up to a remainder TZg = 05(max(/i 2 , /it 2 )). Proposition 
now proved for 9 6 K and A = 0. 

Near identity change of variable. Let us consider 

Cs,a = Ce - nvXdlQ , 

where we recall v = e g^^) » an< ^ w ith satisfying (3.7) (with zero on the right hand side). Then 
Ce,A satisfies 

dtC + d x C + eaiC^C + e 2 a 2 C 2 d x ( + e 3 a 3 ( 3 d x C + - - Ai/)figC - + ^)d 2 x d t ( 
(3.8) +Mea a ((«! - 9vax - \vai)C,dl( + («a - 6^ai)(d x C) 2 ) = 5 (max( M 2 , M e 2 )), 



Again, it is now straightforward though technical to check that, denoting £e,A the solution of (3.8) 
and defining v e x as a function of C,e t \ through (3.5 ), then ((e t \,v ex ) satisfies (3.6 1 up to a remainder 
7lo,X — Cs(max(/i 2 , ^e 2 )). Proposition 



3.1 



is now proved for any 9, A G 



□ 



Discussion. We note that the accuracy of the unidirectional approximate solution, described in 
Proposition |3 . 1 1 is tremendously improved, when compared with the results of section [2.2| Indeed, 



using the stability Hypothesis |2.12| one would obtain the following result. 

Corollary 3.5 (Convergence of unidirectional approximation). For (e, fx, 5, 7) = p G V, as de- 



fined in (1.2), let Uq N be a solution of Green-Naghdi equations (JTTTj) such that the family (Uq N ) 



is uniformly bounded on H s , s sufficiently big, over t ime interval [0, Tqn]> an d with initial data 
satisfying (3.2). Assume that hypotheses of Proposition 3.1 hold, and denote Uq L the unidirectional 



approximation defined therein. Then if Hypothesis 2.12 is true, one has for any t < min(TGw, T s+ §), 



Ucl~Us 



i»([0,i];H») 



< C max(e 4 , pL 2 ) t, 



with C = C( sup„ 



\U P I 
\ u gn\ 



s+5, 1 



c|A|,|^|). 



Such a result is supported by numerical simulations. In Figures [6] and [7J we compute the 
decoupled Constantin-Lannes approximation of Definition |2.3| as well as the unidirectional approx- 
imation described in Proposition |3.1| and compare them with the the solution of the Green-Naghdi 
system (1.1), making sure that initial data satisfies (3.2), and in the Camassa-Holm regime e 2 = \x 



(see Section 2.3 for a description of our numerical method). The display of the results is similar to 



the ones in section 2.3 except we add a plot of the deformation at final time in the bottom panel. 



We also choose to limit the time domain of our numerical simulations to 1/e (instead of 1/e 3 ' pre- 
viously) as more interest is given to short time scales, where decoupled and unidirectional models 
are significantly different, rather than very long time scales, where nothing relevant happens. 



In Figure [6j we choose a set up which is favorable to the CL approximation of Definition |2.3| 
critical ratio 5 2 = 7 and localized initial data. We see that the unidirectional model is still a much 
better approximation. In particular, the CL decoupled model predicts a small wave moving on the 
left (of size C(e 2 ), as ct\ = in (3.2)), which is not predicted in the unidirectional model, and 
almost nonexistent in the solution of the Green-Naghdi system, as we can see in Figure 6(c) This 
short-time 0(e 2 ) error of the CL decomposition is preserved over times of order T = 0(l/e). As 
for the unidirectional model, the produced error is clearly, and as expected, of size C(e 4 t). In the 
Camassa-Holm regime, such accuracy is of the same order of magnitude as the full Euler system 
itself. This means that if the initial data satisfies 



approximate solution described in Proposition 3.1 



3.2), then the very simple unidirectional model 
is as precise an approximation of the solution 



of the full Euler system (A.4) as the solution of the coupled Green-Naghdi system (1.1) (provided 



that the stability Hypothesis 2.12 holds); see Remark 2.14 

In Fig ure [7| the ratio is non critical 5 7^ 7. The accuracy of the CL approximate solution of 
Definition |2.3| is worse than in the critical case, as the short-time error is of size 0(e). Once again, 
the same error estimates holds over times of order T = C(l/e). The accuracy of the unidirectional 
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(a) behavior of the error with respect to time (b) behavior of the error with respect to e = ^fji 
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(c) Error at time t = 20 for e = 0.05 

Figure 6: Unidirectional, localized initial data, critical ratio, Camassa-Holm regime 
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(a) behavior of the error with respect to time (b) behavior of the error with respect to e = -^/ji 
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(c) Error at time t = 20 for e = 0.05 

Figure 7: Unidirectional, localized initial data, non-critical ratio, Camassa-Holm regime 
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model is not affected, and is still of size C(e 4 t): the criticality of the depth-ratio do not play a role 
in the accuracy of the unidirectional approximation. 

Let us now turn to the following question: is it true that after a certain time, any perturbation 



will decompose into two waves, each one satisfying (approximatively) an equation of the form (3.2 I ? 



Our answer is numerical. We use the numerical simulations of section 2.3.2 (non-critical ratio 



localized initial data), and test the right wave of the numerical solution of the Green-Naghdi system 
against (3.2). As we can see in Figure 8(a) (where the log of the error is plotted to ease the viewing), 



a very good agreement appears after a given time T$, which is independent of e (but rather depends 
on the thickness, or wavelength of the initial data). The accuracy of this agreement is in our 
simulation of size 0(e 4 ), and valid for long times; see Figure 8(b) 



epsilon = 0.1 0.05 0.035 




e 




0.035 0.05 0.075 0.1 
epsilon 



0.035 0.05 0.075 0.1 
epsilon 



(a) behavior of the error with respect to time 



(b) behavior of the error with respect to e = ^/Jl 



Figure 8: Validity of (3.2) for generic initial data 



A Derivation of the Green-Naghdi system 



This section is dedicated to the construction and justification of the Green-Naghdi model (1.1 1, 
which is the groundwork of our study. We first recall the so-called full Euler system, governing 
the behavior of two layers of immiscible, homogeneous, ideal, incompressible fluids only under the 
influence of gravity. The derivation of such equation is not new, and we refer to [5] for more details. 

This system can be written as two evolution equations using Zakharov's canonical variables [53] , 
namely £ the deformation at the interface, and the trace of a velocity potential at the interface. 
Such a formulation relies on so-called Dirichlet-to-Neumann operators, solving Laplace's equation 
on the two domains of fluid, with suitable Neumann or Dirichlet boundary conditions. 

We non-dimensionalize the system in order to put forward the relevant dimensionless parameters 
of the system, and in particular e, the nonlinearity parameter, and /i, the shallowness parameter. 

The key ingredient in the construction of the Green-Naghdi model is to obtain an expansion of 
the Dirichlet-to-Neumann operators, with respect to the shallowness parameter, /i. Our result is 
disclosed in Proposition |A.2[ below. 

When replacing the Dirichlet-to-Neumann operators by the truncated expansion, one obtains the 



Green-Naghdi models: (A. 14), (A. 17) and (A. 19) (each version handles different velocity variables 



as unknown. The latter is the system we base our study on, and considers the shear layer-mean 
velocity, obtained after integrating the velocity potential across the vertical layer in each fluid [^] ) . 



2 Note that one cannot write the full Euler system in the simple form of two evolution equations using layer-mean 
velocity variables, as the pressure cannot be eliminated from the equation. The nice formulation of the Green- 
Naghdi system with the shear layer-mean velocity relies on the assumption of shallow water, /i <SC 1, which allows to 
approximate Zakharov's canonical variables in terms of layer-mean velocity variables. See also 1151 and |16| for the 
formal construction of Green-Naghdi models using layer-mean velocity variables from the beginning. 
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All of these asymptotic models are justified by a consistency result, stating that any solution of the 
full Euler system will satisfy the Green-Naghdi systems up to a small remainder (of size 0(p 2 )). 

A.l The Full Euler system 



z 



^-3^ = 



vi = Vt,z0i div Vi = A<f>i = 
f{Pl = 



|2 _ _ P_ 
Pi 




d,Q = x/1 + I d*i ?8n<t>i = V 1 + I 9 *C ?3n<t>2 



v 2 =V ;c , z 02 divv 2 = A0 2 = 9t02 + I |V a; , z (/> 2 \ 2 = -f 2 

3^02 = 



-d 2 



X 



Figure 9: Sketch of the domain 

The system we study consists in two layers of immiscible, homogeneous, ideal, incompressible 
fluids only under the influence of gravity (see Figure[9]). We restrict ourselves to the two-dimensional 
case, i.e. to horizontal dimension d = 1. 

We assume that the interface is given as the graph of a function £(t, x) which expresses the 
deviation from its rest position {(x,z),z — 0} at the spatial coordinate x and at time t. The 
bottom and surface are assumed to be flat. Therefore, at each time t > 0, the domains of the upper 
and lower fluid (denoted, respectively, Vl\ and f2|), are given by 



fil = { (x, z) e 



9J 2 



{ (x,z) e R d x 



C(t,x) < z < di }, 
-da < z < {(t,x) }. 



We assume that the two domains are strictly connected, that is 

di + C(t, x) > h > 0, d 2 + C(t, x) > h > 0. 

We denote by (pi,Vi) and (p 2 , V2) the mass density and velocity fields of, respectively, the upper 
and the lower fluid. The two fluids are assumed to be homogeneous and incompressible, so that the 
mass densities pi, P2 are constant, and the velocity fields Vi, v 2 are divergence free. As we assume 
the flows to be irrotational, one can express the velocity field as gradients of a potential: Vj = V4>i 
(i = 1,2), and the velocity potentials satisfy Laplace's equation 



di 



d 2 z & = 0. 



The fluids being ideal, they satisfy the Euler equations; the momentum equations can be inte- 
grated, which yields the Bernoulli equation in terms of potentials: 



d t <t>i + dV*,. 



P 
P 



--gz inO^ (i=l,2), 



where P denotes the pressure inside the fluids. 

From the assumption that no fluid particle crosses the surface, the bottom or the interface, one 
deduces kinematic boundary conditions, and the set of equations is closed by the continuity of the 
pressure at the interface, assuming that there is no surface tension. [[] 

3 The surface tension effects should be included for our system to be well-posed. However, the surface tension is 
very small in practice, and does not play any role in our asymptotic analysis. See |4.'il for in depth study of this 
phenomenon. 
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Altogether, the governing equations of our problem are the following: 



(A.l) 



d%<j>i 
d t <f>i - 



= 

= 



_p 

Pi 



dtc = 

d,<h = o 

P continous 



\d x C\ 2 d>. 



n<P2 



in n\, i = 1,2, 

in i = 1,2, 

on r 4 ee {(x, z),z = di}, 

on r ee {(x, z), z = £(i, x)}, 

on r b ee {(x, z),z = -d 2 ], 

on r, 



where n is the unit upward normal vector at the interface. 



Rewriting the system as evolution equations. The key point is to remark that the sys- 
tem ( A.l ) can be reduced into two evolution equations coupling Zakharov's canonical variables [63] . 
namely the deformation of the free interface from its rest position, £, and the trace of the upper 
potential at the interface, ip, defined as follows: 

l/j EE (f>l(t,X,£(t,x)). 

Indeed, 4>i and 4> 2 are uniquely deduced from (£, tp) as the unique solutions of the following Laplace's 
problems 

C ( dl + d\ ) fa = in fii, f ( d 2 x + 31 ) 2 = in 2 , 

(A.2) Wi=ii' on T, and < d n (j) 2 = <9„0i on T t , 

{ d z 4>i =0 on T t , [ d z (j) 2 = on T b . 

More precisely, we define the so-called Dirichlet-Neumann operators. 

Definition A.l (Dirichlet-Neumann operators). Let ( e lU 1 ' 00 ^), and d x ip £ fl" 1 / 2 (M). Then we 
define 

G[C]V = v / iT[W(^i) U= c = (^i)U= c -(5xC)(Mi)U=c. 

F[C]V = 0*(&U=c) = 5*(^2(t,a!,C(*,a?))). 
where <f>\ and (f) 2 are uniquely defined (up to a constant for 4> 2 ) as the solutions in H 2 (M) of (A.2). 



The well-posedness of Laplace's problem (A.2), and therefore the Dirichlet-Neumann operators, 
follow from classical arguments detailed, for example, in [22J, Proposition 2.1]. 



One can then rewrite the conservation of momentum equations in ( A.l ) at the interface, thanks 
to chain rule: 

u , li fl /, , m2 , (G[M + (d x Q(d x (<h\,=Q) 2 
d t (^c)+9C+- 2 \d x (H\ z ^)\ + 2{l + \d x tf) = — ' 

Using the continuity of the pressure at the interface, one deduces from the identities above 

d t (p 2 H[C]i> - pid x yj)+g( P2 - Pl )d x C + \d x {p2\H[CW ~ l\d x ^\ 2 ) = W(C», 
where TV is defined as 

jV(CV) = Pl{Gm + (^C )(g^)) 2 - p 2 (G[C]^ + (d x QH[t]if) 2 



2(i + |a,ci 2 ) 



The kinematic boundary condition at the interface is obvious, and the system (A.l ) is therefore 
rewritten as 

(A3) / ^(^[C]^-Pi9^)+5(P2-pi)^C+5^(p 2 |i?[C]^| 2 -pi|^| 2 ) = Wc,V-), 
which is exactly system (9) in [5]. 

4 The solution of the second Laplace's problem is denned up to a constant, which does not play any role in our 
analysis. 
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Nondimensionalization of the system. Thanks to an appropriate scaling, the two-layer full 
Euler system (A. 3) can be written in dimcnsionless form. The study of the linearized system 
(see [43], for example), which can be solved explicitly, leads to a well-adapted rescaling. 

Let a be the maximum amplitude of the deformation of the interface. We denote by A a 
characteristic horizontal length, say the wavelength of the interface. Then the typical velocity of 
small propagating internal waves (or wave celerity) is given by 



[p2 - Pi)did 2 
pid\ + P\d 2 



Consequently, we introduce the dimensionless variable^ 



the dimensionless unknowns 



z — , x — , t — t, 
di A A 



a aAco 



and the four independent dimensionless parameters 



Pi a _ d\ di 

1 = — , e = -J-, M = T2' ° = T- 
P2 «i * d 2 



With this rescaling, the system (A. 3) becomes (we withdraw the tildes for the sake of readability) 
(A.4) 



dtC ~ -G^ = 0, 
P 



with 



9 t (V'^- 7 a^) + (j + 8)d x ( + e -d x (\H^\ 2 - 7 \d x ^\ 2 ) = ped x M^ 

= (iG^yj + e(d x QH^) 2 - j(^G^ + e(d x Q(d x ^)) 2 
2(l + n\ed x C\ 2 ) 

and the dimensionless Dirichlet-to-Neumann operators defined by 

C'fy = vT + ]i^CF(»n^i) U=ec = -Me(^C)(^i) U= eC + (Mi) U=eC, 
ff'V = d x (<h.\ z=t< ) = (^i)| J=£ (+e(«)(^i)| z=£(l 

where cf>± and <p 2 are the solutions of the rescaled Laplace problems 

( pd 2 + d 2 ) 4>x = in fi a = {(a:, z) e R 2 ,e((x) < z < 1}, 
(A.5) <( a z 1= O onr ( E{z = l}, 

^2 = on T = {z = eC}, 

( A*fi£ + flj ) ^2 =0 in f} 2 = {(*,*) £R 2 ,4<K eC}, 
(A. 6) <( 3„0i = S„0 2 on T 2 , 

0,02=0 onr b = {z = -i}. 

Again, the Dirichlet-Neumann operators are well defined, provided that one has ( 6 W 1 
d x ip € i/ 1 / 2 (M), and the following condition holds: there exists h > such that 

(A.7) h x = 1 - eC > h > 0, /i 2 = v + < > & > 0. 





5 We choose di as the reference vertical length. By doing so, we implicitly assume that the two layers of fluid have 
comparable depth, and therefore that the depth ratio S do not approach zero or infinity. Note also that eg — > as 
P2 -> Pi ■ 



28 Decoupled and unidirectional asymptotic models for the propagation of internal waves 



A. 2 Asymptotic models 



Our aim is now to obtain asymptotic models for the system ( A.4), using smallness of dimensionless 



parameters. In the following, we will consider the case of shallow water, namely 

H < 1. 

The key ingredient is to obtain an expansion of the Dirichlet-Neumann operators, in terms of /x. 
Replacing the operators by the first order terms of these expansions allow to obtain the desired 
asymptotic models, which is justified in the sense of consistency (as precisely explained below). 

As a second step, we will rewrite the equations in terms of the shear layer-mean velocities. 
One benefit of such a choice is that it yields to a much better behavior concerning the linear 
well-posedness, as we discuss at the end of this section. 

Our method has been used by Alvarez-Samaniego and Lannes pQ in the case of the water- 
wave problem (one layer of fluid, with free surface), and lead the authors to a complete rigorous 
justification of the so-called Green-Naghdi equations 28J. In the case of two layers with a free 
surface, a shallow water model (first order) and Boussinesq-type models (in the long wave regime) 
have been derived and justified in the sense of consistency in [5], the analysis below is therefore an 
extension of their work. Similar models as our Green-Naghdi system have been formally obtained 
in [16], as well as in [51] (with the additional assumption of 7 ps 1) and in [20]. Let us also mention 
the work concerning the case of two layers of fluids with an interface and a free surface: Green- 
Naghdi-type models have been derived in [SIS], and justified in the sense of consistency in |22) . 
One could formally recover our models from [JH (44) and (60)] by forcing the surface to be flat 
(a = using notation therein). 

Expansion of the Dirichlet-Neumann operators. 

The main ingredients of the following Proposition are given in [5] , and we extend their result one 
order further. 

Proposition A. 2 (Expansion of the Dirichlet-Neumann operators). Let so > 1/2 and s > so + 1/2. 

Let -0 be such that d x ip G H S+11 / 2 (R), and ( £ H s+9 / 2 (R). Let hi = 1 - e( and h 2 = 1/5 + e( such 



that (A.7) is satisfied. Then 



(A.8) |-G"^V ~ 9,(AiW)| H . < M<?o, 



(A.9) |-G^ - d x (hAi>) ~ »l%(hi S d*i>)\ HS < C u 

f-l O 



(A.10) |if^ + -i^U < MO), 

1 h 2 ,H 

(A.ll) \H^i;+^d x ^- ^d^h^d^^d^-h^d 2 ^)]^ <^C X , 

with Cj = C(i,e max ,^ max , j^-,<S max , \C\ Ha+5/2+2] ,\d x ^\ Hs+7/2+2j ). The estimates are uniform with 
respect to the parameters e € (0, e max ), /1 € (0, ^ max ), 6 € (S min , £ max ) and 7 e (0, 1). 

Proof. As remarked in [5], the operator G^'^ip can be deduced from similar operator in the (one 
layer) water wave case with flat bottom: 

= -g[-e<Z]1>, 

where Q is defined in [TJ section 3], and estimates ( A.8 ),( A~9| follows from Proposition 3.8 therein. 



Estimate (A.lOl is given in [SJ section 2.2.2], and we obtain (A.ll I using the same method 



expanding one order further. Let us detail the strategy. 



The first step consists in rewriting the scaled Laplace problem (A.6) into a variable-coefficient 



boundary- value problem on the flat strip 5:=Ex ( — 1,0) using the diffeomorphism 



s ^n 2 , 

a: (x,z) 1— > <j(x, z) = (x, (1/5 + e()z + e£)- 
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Now, one can check (see [5] or [521 Proposition 2.1]) that <f>2 solves (A. 6 1 if and only if cf>2 = 4>2 ° o 
satisfies 



(A.12) 
with 



V x ,*-Q M [cC]V a .,,^ = in 5 

C«2 U=o = d n cj> 2 |, = _i - 0, 



^<9 z ct -/J.d x a 







1/5 + eCV 1 



1/(5 + eC -(z + lJe^C' 

-(, + ivac l(z + 1)£ ^ c|2 



and where d n (f>2 stands for the upward conormal derivative associated to the elliptic operator in- 
volved: 



d n 4>2 \z=z = 



z92 \z=z 



The asymptotic expansion of the Dirichlet-Neumann operator i/ M ' e is deduced from the identity 

= d x U 2 \ z =o)- 



The second step consists in computing the formal expansion of 02 the solution of (A.12), as 

4>2 = 4> (Q) + + v 2 4> (2) + A^r, 



One can solving (A.12 1 at each order, using the obvious expansion of Q M , as well as the known 
expansion of the operator G^tp. This yields explicit formulas for and the estimate follows 
from adequate control of the residual <p r . 



At first order, one has 



jj^d^ U =0 = jji^d^ \ z= _ x = o, 



so that 



<j> w (x,z) = <f> {0) (x), 
independent of z. At next order, one has (denoting /i2(x) = 1/6 + e((x)) 



-8tt<n = -v. 



l>2 

-(z + l)ed x ( 



-(z + l)ed x ( 
|(z + l)e^C| 2 I V^(°) = -h 2 d 2 x 



h 2 z 



in S 

-0, 



where we used that <p > is independent of z. The above system is a second order ordinary differential 
equation, which is solvable under the condition 



d x (h!d x i>) = -d x (h 2 d x <f>W), 



and whose solution is then 



with <pQ (x) being a function independent of z, to be determined later. Note that since the hori- 
zontal dimension is one jf]and using the fact that the fluids are at rest at infinity, integrating the 
compatibility conditions yields 

d x ^ = 

h-2 



6 In the 2d case, one should introduce the non-local operator of orthogonal projection onto the gradient vector 
fields, as in [8], to ensure that the right hand side is a gradient. 
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Let us turn to the next order. One has 

= -h2d 2 A 1} + (z + l) 2 F(x) in S 

i_dA {l) U=o = G{x) + {d x h 2 ){d x ^ ] ) + Wh 3 ^) £0,^°) U=-i = o, 

with 

F(x) = ^fd x ^ and G(x) = -^(h^d^ . 
Solving the first identity with boundary condition d z <f>^ |z=-i = yields 

^(x,z) = - i - Z ^h^d 2 J^F(x) + ^^h 2 F(x)+^\x), 

with (f> (x) independent of z (and which can be set to zero for simplicity) , and solving the boundary 
condition at z = yields the compatibility condition: 

-h 2 d 2 J^ + ±F(x) = G(x) + (d x h 2 )(dJi ] ) + ^(/u 3 c&/>), 

or, equivalently, 

d x {h 2 d x ^ ] ) = l -d x {h^dl^) - ^d 2 x (ht 3 d 2 x i>). 
Finally, integrating this identity and using the expression of d x (f>^ obtained above, one deduces 

The final step is as follows. Let us define 

<fe,app = (O) + 

where (f>^ and (p^ have been obtained by the above calculations. Note that 

ffapp = a c ( fo,app U=o) = d x 4> { -°\x)+ l ld x ( K -^h 2 2 d 2 <f>W +<f>£\x) 



h 



which is exactly the expansion in (A. 11 1. Therefore, the result follows from an adequate estimate 
on 

U = d x ((f>2 U=0 ) ~ d a ( 02,app \z=0 )■ 

Let us first note that one has straightforwardly 

i f max,/ i maj[j |C|#s+9/2i | 9 x 1p | jjs+n/v)- 

Then, using previous calculations, v = <f> 2 — </>2,app — /i 2 ^ 2 -* satisfies the system 



V s ,»-Q^[eC]V»,»« = A* 2 V^,-h in S 
i VU] t d n fc\ z=0 = V + t , 2 (J) .h| 2=0 $.&L=-i=M a (j) "h 

with h = Q^V^^ 2 ', and y = G^tp + nd x {h x d x i)) + fi 2 g<9 2 0i 3 <9 2 V) , so that (TO) yields 



I^Ijjo — A 1 C( ^ i e maxj Mmaxj I C| i I ^a: V'j s + ll/2 ) ■ 

One can now apply Proposition 3 of [8], after straightforward adjustments, and deduce 

\d x V \ z =0 |jjs ^ M ? £ max: Mmaxj | C| ? | | ^s+ll/2 ) ■ 

The Proposition is proved. □ 
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Let us now plug the expansions of Proposition A. 2 into the full Euler system ( A.4), and withdraw 
C(/i 2 ) terms. One obtains 



d t ( - d^hAiP) - ^d 2 x (hi 3 d 2 x ^) = 0, 



(A.14) < 



hi + jh-2 



«2 ' 



2 \ ri2 ri2 
where we denote Vd x ib = V[hi, h2}d x tp, and with the following operators 

V[hxMV = \{h 2 i d x { } ^V)-h l a d x v), 

Af[hi,h 2 ]V = \((d x (hiV)-e(d x 0Y 2 V ) 2 - l(hid 2 ^y 

Our model is justified by the following consistency result. 

Proposition A. 3. Let U = be a solution of the full Euler system (A.4) such that ( |A.7 ) is 

satisfied, and £ e H^([0, T); H s+11 / 2 ), d x ip G W x ([0, T); tf s + 13 /2) s > s + 1/2, s > 1/2, and 
where W^fO, T); iP +1 (IR)) denotes the space of the functions f(t,x) such that 



l/l 



W 1 ([0,T);_H' s + 1 (R)) 



l/l 



L°°([0,T);ff s + 1 (R)) 



l^/l 



L°°([0,T);ff 5 ( 



< oo. 



TTien [/ satisfies (A.14), «p to a remainder R, bounded by 

< M 2 C, 



R 



with 



C - C (s -i/2> S- 



, ^max, A*max, i "max, || S> || yyl rm j>\.^ys+ll/2j j 



L°°([0,T);H') 



W 1 ([0,T);iT s + 13 / 2 ) 



), uniformly 



with respect to the parameters e € (0, e max ), /x <E (0, /i max ), <5 € (Smim <5max) and 7 € (0, 1). 

Proof. When plugging U into ( A.14[ ), and after some straightforward computation, one can clearly 
estimate the remainder using the expansions of the Dirichlet-Neuman operators in Proposition |A.2| 
The only non-trivial term comes from dt (H»< e <ip - ~jd x ^j , which requires the corresponding expan- 
sion of dtH^^tp. This can be obtained as in Proposition |A.2| using the elliptic problems satisfied 
by the time derivative of the potentials. The expansion follows in the same way, provided that 
C € W l {[Q,T);H s+11 / 2 ),d x ij) € W^QO, T); £P+ 13 / 2 ). See the proof of [22 Proposition 2.12] for 
more details. □ 

In [5], the authors use the shear velocity as for the velocity variable: 

(A.15) v = x ((0a-7&)U=, c ) = H^j-jd^. 

The expansion of ff M ' e in ( A.ll| ) allows to obtain approximately d x ip as a function of v: 

h 2 



(A.16) = — h 2 v + fi 

hi + 7/12 n 



~^ 1 -dJv[h 1 ,h2](-j- 

1 + 7/12 V hi 



7/12 



0{f 



Here and in the following, we use the notation O(-) for estimates as in Proposition A. 2 



Plugging (A.15) into (A.14), and again withdrawing C(^ 2 ) terms, yields 

hih 2 



(A.17) 



8 t v + (j + 6)d x ( + ^(7^ 



hi -lh 2 . |2 
lh~2? 



0, 



/j,ed x (TZ[hi,h2]v) 
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with the following operators: 



h 2 V 



•7/12 



ft 2 v 



System (A. 17 1 is justified as an asymptotic model for the full Euler system ( A. 4 1, by a consistency 
result, with the same precision as (A. 14). 



Proposition A. 4. Let U = (CVO be a solution of the full Euler system (A.4) such that (A.7) 
is satisfied, and , .J_€W 1 ([0, T}; H s+11 ' 2 ), d x ipeW 1 ([0, T}; H s + 13 / 2 ) with s > s Q + l/2, s > 1/2. 
Define v as in (A. 15). Then (C>«) satisfies (A. 17), up to a remainder R, bounded by 



R 



L«([0,r);H«) 



With C — ^( s — 1/2 ' h ' £max ' Mmaxj j ^maxj 1 1 C || yj\ (\0,T];H s + 11 / 2 ) ' 1 1 ^ ^ 1 1 W 1 ([0,T] ;.£f s + 13 / 2 ) ) ' un if orm ^V 

with respect to the parameters e G (0, e max ), (i € (0, /i max ) 7 (5 € (<5 m in, <5 m ax) 7 £ (0,1). 

The proof of Proposition |A.4| is identical to the proof of Proposition |A.3| once one obtains a 
rigorous statement of (A. 16) in VF 1 ([0, T]; H s+1 ) norm, thus we omit it. 

Using layer-mean velocities. 

As we have seen, several different velocity variables are natural when expressing the Green-Naghdi 
equations. In the following, we choose to use, as in [T51 [T5] for example, the shear layer-mean 
velocity, defined by 

(A. 18) v = u 2 — 7U1, 

where U\, u 2 are the layer-mean velocities integrated across the vertical layer in each fluid: 



u\(t,x) 



1 



h\(t,x) JeC(t,x) 



d x 4>i(t, x, z) dz, and U2(t,x) 



1 



h 2 (t,x) J_i 



i((t,x) 



d x 4>2(t, x, z) dz. 



We see two main benefits for such a choice. First, the equation describing the evolution of the 
deformation of the interface is an exact equation, and not an C(/i 2 ) approximation. What is 
more, the system obtained using layer-mean velocities have a nicer behavior as for the linear well- 
posedness. These two facts shall be discussed in more details below. 



When integrating Laplace's equation in ( A.l ) against a test function (p{x, z) = tp(x) on the lower 
domain Q 2 , and using boundary kinematic equations, one has 



= 



^d 2 x + d 2 z )cj> 2 dxdz = - [[ V£,,0 • V£,& dxdz + f 



ip d n <j) 2 - 



= — I dxy/Jidxip / \/J^d x (f> 2 (x, z) dz — I dxG^'^tp, 
JM J-l/S Jr 

where V[i x ,z = (^/ft-dx, d z ) T , so that one deduces that for any x G 

1 



-d x {h 2 u 2 ) 



In the same way, one has 



(p{pid 2 + d 2 )<f)i dx dz 



z<P ■ ^xz<t>i dx dz - / ip d n (f>i + ip d z 
1 Jr Jv t 



= — I dx dxy/fitp / y/JId x <f)i(x, z) dz + dxG^'^ip, 
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so that 



0x(toh) = -G^. 



One deduces from the two above identities, imposing zero boundary conditions at infinity (the 
fluid is at rest at infinity), 



7 — , - j _ hi + 7/12 _ 

h 2 u 2 = —hiui, and v = ; u 2 

hi 



hi + 'fh 2 



-Ui, 



It follows that the first equation in (A.4| becomes 



1 



d f ( = -G^ip = -d x 



hih 2 



hi + 7/12 



Let us retell that this identity is exact, as opposed to the C(/i 2 ) approximations in previous asymp- 
totic models. 

One also deduces an expansion of v in terms of v, using Proposition A. 2 (or, equivalcntly, 
identifying the above identity with the first line of (A. 17)). It follows 

hi + 7/12 



v = v + (1- 



h x h 2 



-Q[hi,h 2 }v + 0(p 2 ). 



System (A. 17) therefore becomes 
(A.19) 

hih 2 



7/12 



0, 



dJv + fiQ[h u h 2 ]v\ + (j + S)d x ( + )2 \v\-) = !><<>. (■£[/,, .1, ..}r) . 



with the following operators: 



Q[hiMV 
TZ[hi,h 2 ]V 



^( M -( fe3a -(iiTT^)) + i^'^i 



h 2 V 
7/12 



7^2 



fca V 



+ 



- l(hid x (— 
d x (h 2 3 d x (-^ 



3/11+7/12 V Al2 V V /li+7/l 2 



'/ii V 1 x{ hi+jh 2 J J 



Again, one obtains the consistency of the full Euler system towards our Green-Naghdi model, after 
several technical but straightforward computations. 



Proposition A. 5. Let U = (OVO be a solution of the full Euler system (A.4) such that (A.7) is 
satisfied, and ( g W 1 ^, T); H s+11 / 2 ), d x ip € W X {[Q, T); H s+13 / 2 ) with s > s + 1/2, s > 1/2. 
Define v as in (A. 18 1 or, equivalently, by 



1 
/' 



G^^p = -d x 



hih 2 
hi + 7/12 



TTien (C>^) satisfies (A.19), up to a remainder R, bounded by 



\R\ 



L°°([0,T);H s ) 



with C C( SQ _ 1/2 ' h ' £max > Mmax: j ^maxj || C| | W 1 ([0,T);H B + 11 / 2 ) ' 1 1 1 1 W 1 ([0,T):iJ s + 13 / 2 ) ) ' un ^f orm ^y 

with respect to the parameters e € (0, e max ), /i <E (0, /i max ), <5 € (£ m in, 6 max ) and 7 € (0, 1). 
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Remark A. 6. Proposition A. 5 is one of the steps leading to a stronger justification of system (A. 19) 

same precision), in 



(or any consistent model with the 
achieved in the water wave case in JT^. As described in 
from 



the sense of convergence, as it has been 
Section 6.3], the procedure would follow 



(Consistency) This is Proposition A. 5 



( Convergence ) One proves that the consistent solutions of the full Euler system remain close to 



the exact solutions of (A. 19) (and in particular, one shall prove that the asymptotic model is 
well-posed). Such a result is assumed in Hypothesis 2.12 and is currently under investigation 
by the author and collaborators. 



• (Existence) One proves that smooth, uniformly bounded family of solutions to the full Euler 
system exists. This is ensured by Theorem 5.8 in JJW, provided that a small surface tension 
is added, and that an additional stability criterion is satisfied. 

See Theorem 6.1 in J^ffi for the application of such procedure on the so-called "shallow-water/ shallow- 



water" asymptotic model (which corresponds to (A. 19), when withdrawing 0(/i) terms). 



Linear dispersion relations. The linearized system from (A.17) is 



d t v + (j + S)d x C = 0, 



Let us look for solutions of the form £ = (^e 1 ^ v — v o e l ( kx <*■<*), This leads to the algebraic 
system 

-iujv +ik{j + 5)(° = 0, 



^±^v° = 0, 



which yields the dispersion relation (with u>v° — k(j + 6)£ ) 



35( 7 + S) ' 



This equation does not have any real solution w(fc) if fik ^ ^' ^ us s y s ^ em (A.17) 

linearly ill-posed. 



is 



The linearized system from (A. 19) is 



dtC+^fsd x v = 0, 



Same calculations as above yield the algebraic system 



-iuiQ 



_,_ ik -0 
7+« 



0. 



(v° + ^P^ s - ) v°)+ l k( 1 + 6)C = 0, 



so that = ^ffs v ° an d the dispersion relation is 



2 I , , 2 1 + I 5 



This equations has always solutions: lo± (k) = ± 



velocity variables (A. 19) is linearly well-posed. 



V 71 + ^ fc2 W 



Thus the system using layer-mean 
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B Proof of Proposition |2.5 



Here, we study the well-posedness as well as the persistence of spatial decay in Sobolev norms, of 
the following equation: 

(B.l) (1 — fx 0d 2 )d t u + ea\ud x u + e 2 a 2 u 2 d x u + e 3 a 3 u 3 d x ii 

+ fj, vd\u + fie d x (nud 2 u + i{d x u) 2 ^ = 0. 

More precisely, we prove the assertions of Proposition |2.5| 



The existence and uniqueness of u £ C°([0, T/e); H^ +1 ) n C x ([0, T/e); if*) such that u satis- 
fies (B.l) and initial condition u | t=0 = u° 6 H s s > s > 3/2 has been obtained in [THl Proposi- 
tion 4] (where the authors used slightly different parameters). 

The proof is based on an iterative scheme, which relies heavily on the following energy estimate: 



1 



n < 



1 /2 

(\u(t,-)\ 2 Hs +pfx\u(t,-)\ 2 Hs+1 ) = E s (u)(t) < C(M, \u°\ H ,+i), 



l + 0-i I IHJ 

which is proved to be valid for t <T e = C( so _ 1 3 ^ 2 , M, | u° \ H s+i ) /e. One proves in the same way 



p 11 so — 3/2 



s -3/2' 



,M, lu° 



so that the estimate (2.17) follows. 



In order to obtain the estimates in the weighted Sobolev norms, we introduce the function 
v n — x n u, where u satisfies (B.l). We will prove that the following estimate holds for any 
n,k G N, and for < t < T e>fi = C ( 5^33/2 ,M,n,k, YT J= o \x n ^u°\ H ^ k+2j+1 ) x min(l/e, by 
induction on n G N. 



(B.2) |^"^5 t n|| ioo([0)t) + |^"a fcM || ioo([ot) 



n 



1 13 it 



3=0 



Proposition |2.5| follows straightforwardly. 

The case n = 1 . One can easily check that v\ = x u satisfies the equation 
(B.3) (1 — \i 0d x )dtVi + ea\V\d x u + e 2 a.2V\ud x ii + e 3 a^viu 2 d x u 



with 



R\u] 



+ [i vd\v x + (ie d x (nvidlu + b{d x vi){d x u)) = R[ 
-2/u. 0d x d t u + 3[ii>d x u + \ie ((k + b)ud 2 u + 2i(d x u) 2 ) . 



When taking the (L —) inner product of (B.3) with A V\, one obtains (using that the operator 
A s is symmetric for the (L 2 — )inner product, and d x is anti-symmetric) 

Id 2 

- — (E s (vi)) + ( A s (eaiVid x ii + e 2 a 2 v 1 ud x u + e 3 a^v^d^) , A s Ui ) 

-lie ( (wifigu + t^iX&u)) , A*d x vi ) = ( A s R[u] , A s Vl ). 



Now, we use that u is uniformly bounded through (2.17) for t < T e . 

• Using Cauchy-Schwarz inequality, and Moser estimates: |/<?| ffs < ^( So -i/2 ^ \f\H s Hip ^ or 
s > s o > 1/2, one has 



( A s (eaii>i<9 x u + e 2 ct2V\ud x u + e 3 a^ViM 2 (?;„«) , A s «i ) 

!„.*>) S C( 



5 H 2 ^^. 



So - 3/2 



M,| U °| ffS+2 )(£ s M)' 
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In the same way, and using the definition (E s (-)^ 2 = | • \ Hs + /3/i| • \ Ha+1 > c(M) I ' lff s+1 ' 

|( K s (nvidlu + i(d x vi){d x u)) , A"d x V\)\ < C\v\\ H „ +x (\vi\ H „\u\ H3+2 + \vi\ ffs+ i \u\ B .+i) 

1 „, 1 



< 



• Finally, one can check 

\(A s R[u] , A s Vl )\ < ^\vt\ HS C( So |ti| gJ+a + |am| g . +1 ) < Mg( go _ 1 3/2 ,^h | g g + 3)^ s fa). 

Altogether, one obtains the following differential inequality, valid for all t <T e : 



i|( E -W) 2 < e c ( 1 



1 



.s - 3/2 



M,|u°| HS+2 )(i? s («i)) + ^C( sp ,A/,|tt Q | g . + 3)£'M 



Gronwall-Bihari's inequality (see [55] for example) allows to deduce: 

E s ( Vl )(t) < E s ( Vl )\ t=0 + HC{ l — ,M,\u°\ HS+3 ) (exp(C( l — , M, \u°\ H . +a )et) - l) . 

e s -3/2 1 IJ v \ s -3/2 1 / 

When restricting t to t < T £ ^ = C( ag * 3 * 2 , M, \u°\ H s+2) x min(l/e, l//i), it follows that for any 
s > s > 3/2, 

1 



Hi 



sq — 3/2 1 '"f 



In order to control the time derivative of Vi, we take the (L 2 — )inner product of (B.3) with 
A 2s ~ 2 (dtVi). Estimating each term as above, one obtains 

(E^idm)) 2 < e|«i| H .-i|ft^l g .-iC( ao _* |u| fl .) + fKA 8 - 1 ^, A^dxdtvx )| 

+ ^eC|5 f ui| HS + | w i| HS +i| u | gs ) + |(A s_1 i?[u], A s_1 9 t wi)| 

< E s -\d tVl ) (C( \^,M, y\ HS+1 )E s { Vl )+^C{ \—,M, \u°\ HS+2 )) . 

\ s -3/2 s -3/2 J 



Estimate (B.2) thus follows for n = 1 and k = 0. 

We know turn to xd x u, for fc G N, fc > 1 . Note that 

xd x u = d x (xu) - kd x ~ 1 u, 

so that 



\xd~.u\ 



< \xu\ 



k\u\ 



One deduces, for < t < T(M) min(l/e, !///), 



ii-„ 



< \xu\ 



so - 3/2 



.M, u u 



In the same way, one has xd x dtu — d x {xdtu) — kd x 1 dtu, so that estimate (B.2) hold for n = 1 
and any k € N. 



TTie case n > 2. Let us assume that (B.2 1 holds for any k S N and n < m — 1 (with m > 2). One 
can easily check that u m satisfies the equation 

(B.4) (1 - \x pdl)d t v m + eaiv m d x u + e 2 a 2 v m ud x u + e 3 a 3 v m u 2 d x u 

+ /i vd\v m + d x (Kv m dlu + i(d x v m ){d x u)) = R m [u], 
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with 

Rm[u] = -2n Pmx m - 1 d 3S d t u - 2^ (3m(m - l)x m - 2 d t u + i^mx™- 1 d 2 x u + 3nvm(m - l)x m - 2 d x u 
+ \ivm{m - l)(m - 2)x m ~ 3 u + nemx m ~ l ({K + t)u<9 2 w + 2i(d x u) 2 ) + fj,em(m - l)x m ~ 2 iud x i 

Taking the (L 2 — )inner product of (B.4) with A 2s (v m ), one obtains as above 

l -j t {E s {v m )) 2 < eC( SQ * 3/2 ,M, \u°\ H , +2 )(E s (v m )) 2 + |(A s i^M, A'« m )|. 

Now, one can easily check 

|( A s R m [u] , A s v m )| < fiC{ s j i ,m)\v m \ Hs C(jx m ' 1 d x d t u\ Hs + |: 



^o-l/S 

+ \x m - l dlu\ u , + la^fiLuL. + |x TO - 3 tt 



x m - 2 d t u\ 



m— 1 



Iff' 



< fxE s (v m )c( - l— ,M,m, J2 l gro " 1 " J '»°| g 



As above, Gronwall-Bihari's inequality allows to deduce 



m— 1 

£> m )(*) < £>m)|t=o + C ( So _ 3 /2 ' M ' TO ' S ^"''"^"Iff^+O 



J=0 



for < t < T £ , M = C(^,M,m,E7 = " 
3/2. Therefore, 



m-i | xm _i_ iu oi 



^.s+si+i) x min(l/e, l//z) and any s > s > 



(B.5) \x m u\ H , +1 < <?(— — ,M,m,Y,\x m ~ j u°\ H »+v+i) for < t < T 6 , M . 



The similar estimate on the time derivative of v m , is obtained as above by taking the (L 2 —) inner 
product of (B.4) with A 2s_2 (<9 t w m ). Estimate (B.2) follows for any n = m and k = 0. 

We know turn to x m d x u . Note that for any k € N, k > 1, 



fc-i 



so that 



< \x m u\ H s +k+1 + C(k,m) V \x m - k+j diu\ H s+i. 

j=min(0,fc-m) 



Using (B.5), one deduces by induction on /c that 



1 m 



1 LL 



k— 1 | m— fc+j 



j— min(0,fc— m) 



,s -3/2 



I m -fc +i -i u 0| 



fff, 



i=0 



< C( — ,M,m,k,J2 

Vsq-3/2 ^ 



~™»— j' 01 

a; u B+ k+2j+i 

1 ' u 



for < t < T e ^. 



In the same way, one estimates x m d x d t u by induction, using Leibniz rule of differentiation. 
Estimate (B.2) therefore holds for n = m and any k E N. 



By induction, we proved that (B.2) holds for any k € N and n £ N. Estimate (2.18) follows as 
a direct consequence (using the case k = 0), and Proposition 2.5 is proved. 
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C Proof of Proposition |2.7 



Strategy of the proof. Our strategy is the following. Inspired by the calculations of Section 2.1 
we seek an approximate solution of (1.1) under the form 



app 



= (v + (t,x-t)+v-(t,x + t),(i + 6)(v+(t,x-t)-v-(t,x + t))\ + aU c [v ± ] , 



where v± satisfies the Constantin-Lannes equation ( |2. 12 ), and U c contains the leading order coupling 
terms. The parameter tr is assumed to be small, and we want to justify our approximate solutions 
over times of order 0(1/ a). 

Precisely, our aim is to prove that 

i. the coupling term o~U c can be controlled, and grows sublinearly in time; 

ii. the approximate function J7 app solves the coupled equation (1.1), up to a small remainder. 

As we shall see, controlling the secular growth of U c requires to consider separately the short time 
scale, where the coupling effects may be strong, and long time scale, where the coupling between 
two localized waves moving in opposite directions can be controlled. Thus we introduce the long 
time scale r = t/a, and will seek an approximate solution of (1.1 ) as 

U app {t,x) = U app (at,t,x), with 



U, 



'a.p P (T,t,x) = (v + (r,t, x) + v-(t, t, x) , (<y + 5)(v + (r,t,x) - V-(T,t,x)fj + aU c [v±\(T,t,x) 



with obvious misuses of notation. 

Plugging the Ansatz into the coupled Green- Naghdi equation (1.1) yields at first order 

(C.l) d t v + + d x v + = and 

so that v± (r, t, x) = v±(t,x±) = v±(t, x =Ft). 



dtV- — d x v- 



0. 



At next order, and following the calculations of Section [2Tj one obtains the decoupled equations 

(C.2) er(l - [w t d1)d r v± ± taiv±d x v± ± e 2 a 2 (v±) 2 d x v± ± e 3 a 3 (v±fd x v± 

± /j,(v x + v t )dxV± ± ^ed x (n 1 v ± dlv± + K 2 {d x v±) 2 ) = 0, 

where parameters satisfy identities of Definition |2.3| In order to deal with the coupling terms, we 
introduce the following first order correction. 

Definition C.l. We denote U c = -\-u°_, ( , y+S)(u^_ — ut)j where u±[v±)(T,t,x) satisfies initial 
condition it§_ \ t=o = 0, and equation 

(C.3) a(d t +d x )u c + + f(v + , V-)-f l (v+,0) = 0, and a(d t -d x )u c _ + f (v+, u_)-/ r (0, 5_) = 0, 
where f l and f r are defined as in (2.3a)- |2.3b ). 

Remark C.2. The above discussion does not take into account the use of near identity change of 
variable as described in Section \2.1\ In that case, we set the initial data as follows: 

"» j+S I 



(C.4) 



V± T =0 



(l±/xA^) 



Thus the function v±(at,x^ft) , defined by (C.2) and with initial data ( |C.4[ ) ; satisfies precisely 

v±(at,t,x) = v±(at,x^ft) — v±(t,x^ft), 



2.3 



where v± is the solution of (2.12) as defined in Definition 

The proof is now as follows. First, we state that the approximate solution, U app (o~t, t, x), con- 
structed as above, does satisfy ( 1.1 ), up to a small remainder (actually depending on the size of v±, 
and their derivatives). The fact that v± is uniformly bounded follows from the assumptions of 
the proposition. The core of the proof consists in estimating the growth over long times (that is, in 
the variable t), of u± satisfying (C.3). Each of these steps are described in details in the following. 
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Construction and accuracy of the approximate solution U app . The following Lemma states 
carefully the definition of our approximate solution £/ app , and its precision in the sense of consistency. 



Lemma C.3. Let C°,v° € H s+6 , with s > s > 3/2, and (e, fi, 6,~/) = p G V, as defined in ( |1.2[ >. 
Let v±(r,t,x) be defined by (C.l), (C.2) and with initial data (C.4). Then define v±(r,t,x) as 



v±(T,t,x) = {l±n\dl)- 1 v±{r ) t,x) , 



and set 



U a 



,(r,t,x) = (v+(T,t,x) +v-(r,t,x) , (7 + 5)(v + (r,t,x) - V-(r,t,x))^ + aU c [v±](T,t,x), 



where U c is defined in Definition C.l Then for e small enough, U app {at 1 t,x) satisfies the coupled 
equations (1.1), up to a remainder, 1Z, bounded by 

\fc\ H s < F(e 4 (\v±\ Hs+1 + a\u±\ Hs+1 ) + (J£ 2 {\v ± \ Hs+3 + \d t v± \ Hs+2 + a\u c ± \ Hs+3 + a\u± \ Hs+2 ) 

+o- 2 \d T u± \ R , + eaia\u c ± \ ffs+1 (\v± \ H3 +i+^\u± \ H , +1 ) + e 2 c(| u ± | HS+1 |w± \ HS +i) 
+ 110 {\dlu c ± \ H , + \d 2 x d t u c ± \ H .) (1 + e (\d 3 x v ± \ HS + \d 2 x d t v ± \ H3 )) 

(C.5) + max(e 4 ,^e 2 ,^ 2 )C(|v±| ffs+5 + \d t v±\ HB+i )j , 

with a function F satisfying F(X) < C( sg } 3/2 , M* +6 , 5^-, <5 max , e max , /x max , |A|, \0\)X. 

Proof. This result a direct application of the definitions above, and calculations presented in Sec- 
tion 2.1 More precisely, one obtains 

K = TZ a + <r 2 d T u c ± + (f(v + + au%,v_+au c _)- /(«+,«_)) + Tig +K X , 

where 

• Tvlo is the contribution due to the expansion of h ^]!l / 2 h2 , fej^fe^ ; Q[hi, h,2\v and H[hi, h>i\v 
with respect to parameter e; 



• / is a linear combination of /, and fp^ as defined in ( 2.3a )-( 2.3b); 

• TZg and TZ\ are respectively the components due to the use of the BBM trick and near-identity 
change of variable. 



Let us detail the first terms leading to TZq. We want to control the contribution of the expansion 

hih2 



of u hlh A 1 that is, more precisely, estimate 



TZ, 



(i) 



1 



5 2 - 1 

h 1 +~fh 2 1 + 5 e (i + 5) 2 



:1 5{6 + l) : 
(7 + S) 3 



C +e 



3 7 ^ + l) 2 (l- 7 ) 
( 7 + 5)4 



#3+1 



where we denote C = Ca PP = (v+ + «_ + cr(u+ + u^)(at, t, x), and /ii = 1 - eCa PP , ^-2=5+ eCa PP - 

Note that, as we shall prove that £ app is bounded in L°°([0, T*); H s+1 ), there exists eo > such 
that for any < e < eo, one has 



IClroo <min(l,l/<5) and 



hi + yhi 



< 



Now, one can easily check that, in that case, TZ t 



(i) 



c 4 P {1> «) 
h 1 +~fh 2 ' 



2(1 + 7/5)- 
where pW(X) is a polynomial of 



K 



(i) 



degree 4, and estimate 

As j^rj^ gets multiplied by v and differentiated once, the first term in (C.5) follows. 

/ 2 — h 2 

The contributions due to the expansion of nt^+yfSp ^ s estimated in the same way. The contri- 
bution due to the expansion of TZ[h±, h,2]v requires more derivatives, but may be estimated as above 



40 Decoupled and unidirectional asymptotic models for the propagation of internal waves 



by M e 2 F(| 

v ±\ H s+3 + a \ u ±\ H s+3)- The contribution due to the expansion of Q[h 1 ,h 2 ]v involves 

one time derivative, and is controlled by fie 2 F^jdtv± \ Hs+2 + °i u ± lff s + 2 ) ' This yields the first line 
of < [CT5 ). 

Then, one can check that 

ff»(v + + <ru%,V-+tnf_)- f^(v +) v.) = d x (P(au c ± ,v ± ))+^d 2 d t u c ± 

+ fied x (Qi(au±,d%.v±) + Q 2 {(rd x u±, d x v±) + Qz{ad 2 x u c ± ,v ± ) + Qi(au c ± , d 2 x u c ± ) + Q 5 (ad x u c ± ,d x v ± 

+ ycd t \Q6{(TU±,dlv±) + Q 7 (ad x u c ±1 d x v ± ) + Q s (adlu c ±1 v±) + Q 9 (au c ±1 dlu c ± ) + Q w (ad x u± , d x v± 

where P is a bivariate polynomial of degree 4, and whose leading order terms are 
P(cru±,v±) — aie(au±v±) + e{<ju c ± ) 2 + a2e 2 au±v± 2 + 



and Qi are bilinear forms. Each of these terms are bounded as in (C.5). 

Now, IZg and 1Z\ depend uniquely on decoupled terms v±, and one has similarly 

\fce\ H s + \R-\\ H a < Comax(e 4 ,^e 2 ,^ 2 )(|w ± | ffs+5 + \d t v±\ Ha+i ), 

with C = C( So _ 1 3/2 , j^r, S max , e max , ^ max , |A|, \9\). The proposition follows. □ 



Our aim is now to estimate each of the terms in (C.5 1. 



Estimates of the decoupled approximation v±. The bound on v± is a direct consequence of 
the assumptions of the Proposition, using the identity 

v±(ot, t, x) = v±(at,x^ft) = v±(t 1 x^ft), 



where v±(t, x^ft) is defined through Definition 2.3 and therefore is uniformly controlled as assumed 
in the Proposition. It follows 

(C.6) \\ V ±\\L°°(lO,aT*)x[0,T*);H*+ e ) + HI drV± II L°° ([0,trT*) x [0,T*);H'+S) - M *+6' 

However, let us note that one can gain extra smallness on ad T v± (trading with a loss of deriva- 
tives) from the fact that v± satisfies (C.2). Indeed, one has for any k £ N, 

<r\d T v±\ HS+k < max(eai,e 2 ,^)P(|w ± | H3+fc+3 ) + HVt<r\<%d T v±\ H , +h , 

where P(X) is a polynomial, so that Q 

(C7) cr ll a r S ±IL«([Q,<rT*)x[o,T*);H'+3) - niax(eai, e 2 , /i)C(M* +6 ). 

As for the case of localized initial data, the assumption of the Proposition yields 

6 5 

(C.8) ^\\(l+x 2 ) a d*v4 L „ ([0>aT$)x[m ^ 

k=0 k=0 

and one deduces as above, for k = 0, 1, 2, 3 

a\ (1 + x 2 ) a d k x d T v ± \ H . < max(eai, e 2 , fi)P fe K 1 + z 2 ) Q ^ +i ?) ± \ H )j + H( J + x 2 ) a d k x +2 d T v ± \ H , 
(C.9) < max(ea 1 ,e 2 , A/ )C(Mf +6 ). 



7 Here and in the following, we do not explicitly keep track of the dependence with respect to all the parameters; 
as (e, ft, <5,7) = p £ V , as defined in |L2t , and parameters satisfy identities of Definition |2.3| one should replace 

C(M* + g) by M* +g , <5 ma x, emax, Mmax, |A|, \6\). 
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Control of the secular growth of the coupling terms Let us now study the term 

U c = (u c + + u°_, (7 + 8)(u c + — u c _)) ) which contains all the coupling effects between the different 
components. Our aim is to control the secular growth of this term. This will be achieved through 
the following two Lemmas. 

Lemma C.4. Let s > 0, and f G H S+1 (R). Then there exists a unique global strong solution, 
u(t, x) € C°(R; H s+1 ) n C^R; H s ), of 

(d t + Cl d x )u = d x f with { ( d t + c ^)f = ° 
u\t=o = " " ' \ fi \ t =o = f° 

where c\ 7^2. 

Moreover, one has the following estimates for any t G R: 

2 

')|ffo+l(R) - | Cl _ ~\ \ f lffH-i(R)> l U (*> ')lffs(R) - 1*1 \f Ih«+1(R)' 

Lemma C.5. Let s > sq > 1/2, and v®, G H S+1 (R). Then there exists a unique global strong 
solution, u G C°{R;H S+1 ) n C^ffi;^), of 

(dt+cd x )u = g( Vl , V2 ) v . e{12} Ud t +c i 8 x )v i =0 

u \ t =o =0 { v% | t =o = u" 

where c\ ^02, and g is a bilinear mapping defined on R 2 and with values in R. 
Moreover, one has the following estimates: 



IML~([0,t);ifs) - *l| V llH»(M)r2|fl-»(K)) ll^* u IL~([o,t);H 3 ) - I* I \ V ° I H -+ 1 (K) I W 2 | ff a+i ( 

ii. Using c\ ^ C2, one has the sublinear growth 



^ L.oi L.OI |«„./ + M _^ l + I 01 



01 



V l C l — C 2| V l Cl _ C2 I 

witt c so = c( so _ x 1/2 ). 

m. J/ moreover, there exists a > 1/2 sitc/i t/ia£ + £ 2 ) Q ; ^2(1 + a^ 2 ) a G i/ s (R), then one 
has the (uniform in time) estimate 

\\ u \\l°°hs(r) - ^0 |^i(i + x2 ) a \ H -(m\ V2<yl + a;2 ) Q lff s (R)' 

witfi C So = C'( so j l 1/2 , 0-1/2)- 
Proof. Lemma |C.4| is straightforward, since one has an explicit expression for the solution: 

u{t,x) = -J—(f(t,x+( C 2^ Cl )t)-f(t,x)) = —*—(f{x- Cl t)-f{x-c 2 t)). 
c 2 - Ci c 2 - Ci 

It follows \u\ „,,, < 1 — - — r f° „,.,, and 

I \H 3 + 1 — I C2 — ci I I J \H"+ 1 ' 

' ' ' aj^-y)^ < r \d x f°(x-y)\ Ha dy < \t\\d x f\ Ha . 



u 



C2-Cl\'Jc 3 t * \C2~Ci\J C2t 



As for Lemma [C.5[ the well-posedness as well as estimate i. are standard; the remaining esti- 
mates, controlling the secular growth, are proved in [32], Propositions 3.2 and 3.5. □ 



42 Decoupled and unidirectional asymptotic models for the propagation of internal waves 



Lemma C.6. Let v± defined as in Lemma\CJ^ for (r,t,x) £ [0,crT*) x [0,T*) x R. Then there 
exists a unique U c (t, t, x) £ L°°([0, T) x M; H s ) such that U c = {u c + + u'L, (7 + — «-)) > with 

u± satisfying (C.3) and U c \t=Q = 0. Moreover, one has the estimate for and t £ [0, T*] 



(CIO) 
(C.ll) 



a\\u 



±L»([o,<7T*]x[o,t];H«) - C(AC+3) max ( e «i: e2 ^) mia(i, Vf), 
^ll^rMi Lco ([ o ) r] >< [o )t ];Ha ) < C(M* +6 ) max(eai, e 2 , min(t, V~t). 



Moreover, if (1 + x 2 )(°, (1 + x 2 )u° € H s+3 , s > s > 1/2, and C°,v° £ H s+3+2m , then one has 
the uniform estimate 



(C.12) 
(C.13) 



a \\ u ±\\L°°ilO,aT»MQ,Tt];H») - ^(Mf +3 ) max(eai , e 2 , ^) min(<, 1) , 

CT l^ U ±IL~([0,T]x[0,T»];ffs) - C ( M !+e) max ( m l . <? min (^ X ) 



Proof. The function are defined as the solutions of a transport equation, with a known and 
controlled forcing term, thus the existence and uniqueness of u± and U c is straightforward. 

As for the estimates, we focus below on u c + . Estimates for u c _ follow in the same way. 

By definition, it^ satisfies 

a(d t + d x )u c + = f(v+,0) - f(v+,v-), 
and the right-hand-side can be decomposed in the following way: 

3 3 j 

f(v+,0) - f l (v + ,v.) = y^a.rUAr 1 -) ■ ^«^"-^(''' ri "'i ■ 

3=1 3=1 i=l 

3 

+ ^d x ( Cl v^d 2 x v^ +c 2 (d x v-) 2 ) + ^Y,d k {d k x v+){dl' k v.) 



k=0 



3 3 j 

j=l j=l i=l 



ttd x f 4 + ned x f 5 + A^y^ 



9k, 



fc=0 



where the parameters dj, dij, bj, Cj and dk depend on 7, S, A, 9, and where we used dtV- = d x i 
Note that the first order terms ea± and ea\ : \ are factored by a± — - s ~ 7 



2 7+<5 ' 



The contribution from /.,■, j = 1, . . . , 5 will be estimated thanks to Lemma C.4 as we will use 
Lemma C.5 for the contribution of gij and gk- 

For all < j < 4, one has dtfj — d x fj = 0, and 

• \fj\ H s+i < C\v-\ Lao \v_\ Bs+1 < C(\v-(t,-)\ h . +1 ), for j = 1,2,3; 

• |/4| HS+1 < C\dlv_\ Ha+1 < C(\v-{t,-)\ h , +3 ); 

• l/sl^+i < C|^| L oo|9 2 W-| ffa+1 + |(^v-) 2 | ff3+1 < C(|i5_(T,-)| H .+a). 

It is then straightforward to check that gij and gk are the sum of bilinear mappings, applied to 
functions u±, satisfying dtu± ± d x u± = 0, with u± bounded as below: 

. \vi\ HS < C(\v ± \ HS ) < C(\v ± (t,-)\ h .), fori = 1,2,3. 

• \ 9 x v ±\ H * - c i\ v ±\ H s+i) < C(\v±(t,-)\ Hs+1 ), for j = 1,2,3. 

• |^±| HS < fxC(\v ± \ Hs+k ) < C(\v±(t,-)\ Hs+3 ), for k = 0,1,2,3. 



It follows form Lemmata C.4 and C.5 that u'L satisfies 



o-\u c + {T,t,-)\ HS < min(t, Vt)max(aie, e 2 ,//,/ie)C(|w±(T, Ol^.+a) 
and the last term is uniformly estimated through (C.6), so that ( C.10[ ) follows. 
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As for the case of localized initial data, we make use of the fact that for k = 0, 1, 2, 3, one has 
\(l + x 2 rd k x v ± (r,t,-)\ Hs = \(l + x 2 Td k x v ± (T,0,-)\ Hs = \(l + x 2 rd*vi(r/*,-)\ HS < Mf +3 , 



so that (C.12) follows in the same way for r < <jT^ . 



The estimates on d T v±(r, x) are obtained similarly, as 

3 3 j 3 

(d t + d x )d T u c + = e j d x d T fj + ^ e> ^ d T g id + ^d x d T f 4 + ^ed x d T f 5 + ^e^d T g k . 

3=1 3=1 i=l fc=0 



One can check that each term of the right hand side satisfies the hypotheses of Lemmata C.4 or C.5 



and estimate (C.7) allows to obtain (C.lll 



The case of weighted Sobolev spaces (C.13) follows as above, using estimate (C.9). The Lemma 
is proved. □ 



Completion of the proof. The consistency result stated in our Proposition is now a straight- 



forward consequence of Lemma C.3 together with estimates (C.6)-(C.8|, and Lemma C.6 



One can check that the remainder in can be estimated as 

\K\ H , < C( A C+6)( max(a 2 e 2 ,e 4 ,/i 2 )min(i, Vi) + max(e 4 , ,u 2 ) ) 
Note that we use 

dt{y±(at 1 t, xj) — ad T v± + dtv±, and a\d T v±\ HS < m&x(eai, e 2 , fie)\v±\ HS 
as above, as well as a uniform estimate 

HI u ± IL-([o,r*;H-+i) ^ M - 



The latter estimate can be enforced using Lemma |C.6| by restricting the time interval with 



M 



max(eai, e 2 , fx) 



although such condition is not necessary, as Proposition 2.7 is empty for t > M (max(eai, e , /i)) , 



the accuracy of the uncoupled solutions being of order 0(1). 
More precisely, we proved the following Lemma. 



Lemma C.7. For (°,v° e H s+6 , with s > s > 3/2, and (e,(j,,S,~f) = p G A, as defined in 1.2, let 
U app be defined as in Lemma C.3 Then U ap p(o~t,t, x) satisfies the coupled equations (1.1 1, up to a 



remainder, TZ, estimated for t G [0, T*) by 



L°°([0,t) 



. H s) < C( M s+6){ max(a 2 e 2 ,e 4 ,/i 2 )min(i, Vi) + max(e 4 ,/i 2 ) j. 



and on [0, T') by 



L°°([0 



it);ff .) < C(M» +6 ) ( max(a 2 e 2 ,e 4 ,/i 2 )nfin(M) + max(e 4 , M 2 )). 



The Proposition is a consequence of the above Lemma, together with the estimates of LemmafC.6| 



(see also Remark C.2) 
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